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FOREWORD 


Oceanic  volume  reveiberation  can  adversely  affect  the  peilonnance  of  Navy  sonars.  Small  swimbladder- 
bearing fish  are  the  primary  cauae  of  this  reverberation  at  ship  sonar  frequencies.  This  report  describes  a new 
acouslic  rnodel  of  a swirnbladdec-bearing  fish.  This  nwdel  is  an  irnprovernent  over  previous  rnodels  and  should  be 
of  value  in  future  volume  reverberation  studtos. 

This  report  was  originaly  a dbsertaiion  submitted  in  partial  fulfillment  of  the  requirements  for  the  Ph.D.  degree 
in  acoustics  at  the  Cafhoic  Univetsiiy  of  America. 
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EXECUTIVE  SUMMARY 

T’  A new  model  of  a swimbladder-bearing  fish  has  been  devetoped  in  order  to  provide  improved  predictions  of  the 
resonant  frequency  and  acoustic  cross  section  of  such  a fish.  The  model  consists  of  a small  spherical  shell  in 
water,  enclosing  an  air  cavity  which  supports  a surface  tension.  The  shell  is  a viscous,  heat-conducting  Newtonian 
fluid,  with  the  physical  properties  of  fish  flesh.  A comparison  of  the  results  obtained  with  the  new  model  to 
experimental  data  indicates  that  the  new  model  constitutes  a definite  improvement  over  previous  models.  The  new 
model  can  predict  the  high  values  of  damping  and  elevated  resonant  frequencies  that  previous  models  could  not. 
The  model  appears  to  most  accurate  for  fish  in  which  tension  in  the  swimbladder  wall  has  a minor  effect  on 
resonant  scattering.  This  includes  the  fish  which  are  of  interest  in  studies  of  volume  reverberation  and  therefore, 
the  new  model  should  be  of  considerable  value  in  such  studies. 

^ \ 


D D C 

FEB  14  1978 


D 

lit 


LirLi  li  LL 

D 


I 


ACKNOWLEDGMENTS 

It  is  my  pleasure  to  acknowledge  the  guidance  and  advice  provided  by  Dr.  Thomas  J.  Eisler.  My  appreciation 
also  goes  to  Drs.  Ronald  New  arxl  John  J.  McCoy  for  their  critical  review  of  the  manuscript.  Special  thanks  goto  Mr. 
Robert  S.  Winokur,  for  his  cooperation  and  support  throughout  the  course  of  this  work. 


CONTENTS 


I.  INTRODUCTION  

Volume  Reverberation 

Existing  Models 

Experimental  Data 

II.  A NEW  MODEL 

The  Model 

Limits  and  Physical  Properties 

Formulation  of  Equations  

Boundary  Conditions 

III.  SOLUTION  

IV.  DISCUSSION 

Results 

Comparison  to  Free  Bubble 

Comparison  to  Experimental  Data 

V.  CONCLUSIONS 

REFERENCES  

APPENDIX  A;  Physical  Properties  of  Fish  

APPENDIX  B:  Simplified  Expression  for  (B»/A) 


Page 
. . . 1 
....1 
....2 
....4 
....7 
....7 
....7 
....8 
...16 
...19 
...23 
...23 
...25 
...28 
...33 
..  .35 
...37 
...42 


LIST  OF  ILLUSTRATIONS 


Figure  1 . Resonant  frequency  of  a prolate  spheroidal  butjble 5 

Figure  2.  Resonant  frBquerx:y  of  the  new  model  of  a swimbladder-bearing  fish 24 

Figure  3.  Viscous  arxJ  radiation  damping  factors  for  the  new  model  of  a swimbladder-bearing  fish 26 

Figure  4.  Thermal  damping  factors  for  the  new  model  of  a swimbladder-bearing  fish 27 


LIST  OF  TABLES 

Table  I.  Physical  Properties 8 

Table  II.  Experimental  Data  and  Results  of  Comparisons  to  the  New  Model 29 


IV 


CHAPTER  I 


INTRODUCTION 


Volume  Reverberation 


When  sound  is  propagated  in  the  ocean  any  inhomogeneities  in  the  medium  will  scatter  a portion  of  the 
sound  incident  upon  them.  This  scattered  sound  is  termed  volume  reverberation.  If  volume  reverberation  levels 
are  high,  the  operation  of  an  active  sonar  can  be  adversely  affected.  Specifically,  high  reverberation  levels  can 
mask  an  echo  reflected  from  a particular  target  of  interest. 

During  World  War  II,  researchers  studying  sonar  echo  ranging  at  the  University  of  California  Division  of  War 
Research  (UCDWR)  discovered  that  volume  reverberation  levels  were  usually  vertically  stratified.  Layers  of 
high  reverberation,  on  the  order  of  1 00  m thick,  were  found,  usually  within  the  upper  1 ,000  m of  the  water  column, 
extending  over  large  geographic  areas  of  the  Pacific  Ocean  [1 ,2].  Other  researchers  subsequently  found  such 
layers  in  other  oceans  (3,4).  These  layers  came  to  be  known  as  deep  scattering  layers  (DSL).  In  addition,  it  was 
found  that  these  layers  frequently  rose  to  shallower  depths  around  sunset  and  descended  around  sunrise.  This 
gave  rise  to  the  hypothesis  that  deep  scattering  layers  were  caused  by  biological  organisms,  which  were  known 
to  undergo  diurnal  vertical  migrations  [5]. 

The  great  difference  between  the  acoustic  impedances  of  air  and  sea  water  causes  an  air  bubble  to  be  a 
much  more  effective  scatterer  than  other  objects  of  comparable  size  [6].  This  fact  led  Marshall  to  examine  the 
possibility  that  small  mid-water  fish  which  contained  air-filled  swimbladders  were  the  cause  of  DSL  [7] . His  study 
strongly  implicated  such  fish  as  major  components  of  DSL.  The  primary  scattering  mechanism  was  considered 
to  be  swimbladders  which  could  resonate  in  the  furxlamentai,  or  volume  pulsation,  nxxlevifhen  insonified  at  the 
proper  frequency.  Subsequent  research  by  Mersey  and  co-workers  displayed  the  frequency  dependence  of 
reverberation  levels  and  provided  further  qualitative  proof  that  resonant  scattering  by  swimbladder-bearing  fish 
was  the  major  cause  of  volume  reverberation  in  the  ocean  [8-11].  The  study  of  volume  reverberation  in  the  world 
ocean  has  continued  (for  example,  see  reference  1 2)  and  it  is  now  generally  accepted  that  swimbladders  of  fish 
are  the  predominant  scattering  mechanism  in  most  geographic  areas.  However,  it  has  not  been  until  very 
recently  that  quantitative  comparisons  between  fish  distribution  data  and  acoustic  volume  reverberation  data 
have  been  made. 

The  fundamental  volume  reverberation  parameter  is  the  back-scattering  coefficient  of  a unit  volume  of 
ocean,  M,  which  is  the  ratio  of  the  scattered  intensity  at  a unit  distance  from  the  unit  volume,  I,,  to  the 
incident  intensity,  Ij,  [13], 

M = l,/l,  . (1-1) 

If  it  is  assumed  that  the  scattered  signals  add  incoherently  [13],  then  M can  be  defined  In  terms  of  the  scatterers 
as  n 

I-' 

where  n is  the  number  of  scatterers  in  volume  V and  a,  is  the  ratio  of  scattered  power  to  incident  intensity  of  the 
jth  scatterer.  o is  called  the  acoustic  cross  section  of  the  scatterer. 

Although  M is  the  fundamental  parameter,  the  quantity  which  is  most  often  utilized  in  any  discussion  of 
volume  reverberation  is  the  scattering  strength  per  unit  volume,  S,,  which  is  M expressed  in  decibels. 

S,  = 10  log  M . (1-3) 

Sv  is  often  simply  called  scattering  strength  [14]. 

In  any  acoustic  measurement  of  volume  reverberation,  S.  is  obtained  basically  from  the  equation 

S,  = 10  log  (1,/U  . (1-4) 

The  determination  of  I,  and  I,  is  relatively  straightfonrvard  and  acoustic  measurements  of  Sv  are  fairly  routine. 
However,  If  a determination  of  S„  is  to  be  made  from  the  distribution  of  scatterers,  the  equation  employed  is 
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The  difficulty  in  making  quantitative  comparisons  between  acoustic  measurements  of  volume  reverberation  arxj 
the  distnbution  of  swimbladder  fishes  lies  in  the  determination  of  both  n and  o|.  The  determinaton  of  the  number 
of  fish  arxj  the  species,  size  arxj  swimbladder  size  of  each  is  strictly  a biological  problem  which  will  rx)!  be 
considered  here.  However,  the  difficulty  of  this  problem  is  rx)t  to  be  minimized.  Given  the  size  arxj  swimbladder 
size  of  an  irxjividual  fish,  calculation  of  o is  also  difficult,  due  to  the  complex  structure  of  fish. 

Simplified  models  have  been  developed  in  order  to  estimate  o near  resonarxe  for  an  individual  fish.  These 
models  are  based  on  the  premise  that  only  the  swimbladder  is  a significant  contributor  to  o near  resonance. 
Experimental  eviderx»  [15]  shows,  in  fact,  that  this  is  so,  whereas,  at  frequerxxes  much  higher  than  the 
resonant  frequerx:y,  the  swimbladder  and  the  body  of  the  fish  contribute  about  equally  to  o (16). 

Experimental  eviderx^e  indicates  that  the  existing  models  have  some  shortcomings.  After  a review  of  the 
models  arxj  the  experimental  data,  it  will  be  the  purpose  of  this  report  to  develop  an  improved  model 
of  resonant  scattering  from  an  individual  swimbladder-bearing  fish  in  order  to  eliminate  or  at  least  decrease 
the  si  .vttlcomings  of  the  existing  models. 


Existing  Modeis 

A swimbladder  is  essentially  just  an  air  bubble  within  the  fish,  so  that  the  simplest  acoustic  model  of  a 
swimbladder  fish  is  an  ideal  spherical  air  bubble  having  the  same  volume  as  the  swimbladder.  The 
frequency  of  the  fundamental  mode  of  resonance  of  a small  ideal  spherical  air  bubble  in  water  was 
determined  by  Minnaert  [17]  to  be 

<^23^=2^.  (1-6) 

POw 

where  (lib  is  the  circular  frequency  of  resonance,  a the  equilibrium  bubble  radius,  Pq  the  ambient  pressure,  Va 
the  ratio  of  specific  heats  of  air,  and  Po,.  the  density  of  water.  The  acoustic  cross  section  of  a small  ideal 
spherical  air  bubble  in  water  is  [18] 


o = 


4ntf 


(1-7) 


where  o)  is  the  insonifying  frequency  and  c»  is  the  sound  velocity'  in  water.  The  limiting  factor  on  these 
equations  is  the  size  of  the  bubble,  which  is  limited  to  values  such  that  (u>a/c»)«  1 . 

For  a small  real  air  bubble  in  water,  the  above  equations  for  Ub  and  o remain  the  same,  except  that  the 
term  (cj^a^/c**)  in  equation  1-7  is  replaced  by  d*  [18].  d is  an  unspecified  damping  constant  which  includes 
the  effects  of  heat  conduction,  surface  tension,  viscosity,  and  other  processes. 

Devin  [19]  studied  the  damping  at  resonance  of  real  air  bubbles  in  water.  He  found  that  the  damping 
constant  at  resonance,  6,  was  the  sum  of  three  damping  processes:  thermal  damping,  5»,,  viscous 
damping,  5y„,  and  radiation  damping,  6,^^, 

^ = ^nd  + ^ ■ (1-8) 

6 is  defined  such  that,  at  o)  = (Dq,  d = 5.  Devin  determined  the  values  of  6,*,,  6,^,  and  6,^,  for  bubbles  of  the 
size  which  are  of  present  interest,  to  be 

6r«,=^.  (1-9) 


6 - . 


(1-10) 


where  ri,„  is  the  shear  viscosity  of  water,  x,  is  the  thermal  rxjnductivity  of  air,  pb,  is  the  density  of  air,  and  Cp, 
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is  the  specific  heat  at  constant  pressure  for  air.  Eller  (20)  extended  these  results  to  non-resonant 
frequencies.  However,  Fairbank  [21]  has  recently  contended  that  Eller’s  results  are  incorrect  by  a factor  of 

(t«Jb/a))2. 

A progression  of  improvements  have  been  made  to  Minnaert's  original  work,  which  has  included  such 
effects  as  compressibility,  surface  tension,  and  viscosity  (22-24).  The  resonant  frequency  of  a spherical  air 
bubble  in  water  where  the  effects  of  surface  tension  and  viscosity  are  included  is  (23) 

i.i  2^2  = ^ 2s(3y,  - 1) 

Pow  Pow^3^  Pow® 

where  s is  the  surface  tension.  In  all  of  the  above  references  the  bubble-water  system  was  modeled  as  a 
spring-mass  system  and  the  effect  of  viscosity  included  in  the  damping  term.  Hsieh  (25)  maintains  that  this 
method  of  including  viscous  effects  is  adequate  only  so  long  as  the  damping  is  small. 

In  an  attempt  to  improve  on  the  air  bubble  in  water  as  a model  of  scattering  from  a swimbladder  fish, 
Andreeva  (26)  developed  a model  of  an  air  bubble  in  an  infinite  elastic  medium.  This  elastic  medium  was 
described  by  a complex  shear  modulus,  p,, 

M,  = Mr,  + iM„  (1-13) 

It  has  been  determined  that,  for  fish  flesh  (27) 

10*  dynes/cm^  < m S 10^  dynes/cm^, 
and  0.2  Mr,  <Mrt  <0.3  Mr,, 

Andreeva  chose  Pt,  = 10®  dynes/cm  ^ and  all  other  properties  of  the  elastic  medium  were  set  equal  to  those 
of  water. 

Using  this  model,  Andreeva  obtained  the  following  equations  for  the  resonant  frequency  and  acoustic 
cross  section  at  resonance  of  a small,  spherical  bubble; 

2a2  = 2YaPo_f_^  (1-14) 


/_(Vi_i  ) ' ^ ( 0)0^  \ 

V 0)2  / V ) 


where  Q defines  the  sharpness  of  the  resonance  peak.  Q is  termed  the  damping  factor  and  is  defined  at 
u)  = 0)  0 to  be 

where  o),  and  Wj  are  the  frequencies  at  which  the  scattered  sound  power  is  one-half  that  at  u)o  (19). 
Andreeva  has  utilized  Devin's  (19)  work  to  obtain  values  for  Q. 

1 o)  a 

(I-17) 


_L=  _lik_ 

Q Po,‘*)o^a2  ■ 


1 = 3(y.  - 1)  / K, 

Q ih  a \ 2a)oPo,Cp.  / 


Q was  calculated  by  Andreeva  to  vary  between  5 and  10,  the  value  varying  with  depth.  Equations  1-17  and  I- 
19  are  equivalent  to  equations  1-9  and  1-11. 

Equation  1-14  reduces  to  the  resonant  frequency  of  an  air  bubble  in  water  (equation  1-6)  as  p^ 
approaches  zero.  In  addition,  if  Pp  approaches  zero,  that  is,  if  the  air  were  evacuated  from  the  sphere, 
equation  1-14  reduces  to  the  resonant  frequency  of  a spherical  cavity  in  rubber  (28). 

Equation  1-15  is  actually  defined  only  at  resonance  because  Q is  defined  only  at  resonance.  However, 
for  the  values  of  Q obtained  by  Andreeva,  the  Q-term  is  significant  only  near  resonance.  Therefore, 
although  the  equation  is  not  strictly  defined  away  from  resonance  it  can  be  used  equally  well  at  or  off 
resonance. 
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Lebedeva  [29]  improved  Andreeva's  model  by  considering  an  air-bubble  inside  an  elastic  shell  in  water 
and  by  extending  the  results  to  off-resonance  cases.  However,  thermal  losses  were  neglected.  For  this 


model. 

- T" " ^ 

i'-w] 

1 

(1-20) 

and 

4na^ 

(1-21) 

u)a  , 

('-1^)  ]'■ 

c, 

■p^po^ 

where  b is  the  outside  radius  of  the  shell.  Thus  for  a3/b3<<1,  equation  1-20  is  equivalent  to  equation  1-14 
Equation  1-21  can  be  written  in  the  form  of  equation  1-15  as 

4na^ 
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where,  at  o)  = cjo.  H = Q such  that. 
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In  all  of  the  above  studies  only  the  fundamental,  volume  pulsation,  mode  has  been  considered  and  all 
higher  modes  have  been  neglected.  Strasberg  (30)  has  shown  that  appreciable  sound  is  radiated  from  a 
bubble  only  when  it  is  oscillating  in  the  fundamental  mode  and  the  sound  radiation  associated  with  higher 
modes  is  several  orders  of  magnitude  lower.  On  this  basis,  it  is  evident  the  consideration  of  only  the 
fundamental  rtxDde  is  reasonable,  and  therefore  only  the  fundamental  mode  will  be  considered  in  this 
report. 

In  addition,  all  of  the  studies  have  assumed  a spherical  bubble  or  swimbladder.  However,  swimbladders 
generally  are  not  spherical  in  shape.  They  usually  more  closely  resemble  prolate  spheroids,  with  minor-  to 
major-axis  ratios,  e,  generally  ranging  between  0.1  and  0.5  [6,31-33].  Strasberg  [34]  and  Weston  [35]  have 
shown  that  spheroidal  bubbles  have  a resonant  frequency  only  slightly  higher  than  spherical  bubbles  of  the 
same  volume.  For  a prolate  spheroidal  bubble,  the  resonant  frequency  aio^ , is  [35] 
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(It  should  be  noted  that  this  equation  is  printed  incorrectly  in  reference  35.)  Equation  1-25  is  plotted  in  Figure 
1.  The  assumption  that  the  swimbladder  shape  can  be  modeled  by  a sphere  is  a reasonable  one  for 
1 < E-’  < 4,  where  1 < 2;  <1.1,  and  is  not  really  grossly  in  error  at  an  extreme  e ' of  20,  where  C,  = 1.4. 

Weston  [35]  has  presented  a comprehensive  review  of  scattering,  covering  the  information  available  to 
1966.  Although  a specific  scattering  model  is  not  presented,  many  aspects  of  scattering  from  a fish  with  a 
gas-filled  swimbladder  are  discussed. 


Experimental  Data 

Several  sets  of  experimental  measurements  of  the  acoustic  characteristics  of  swimbladder  fish  near 
resonance  have  been  conducted.  Some  measurements  have  been  conducted  for  the  purpose  of  studying  the 
effect  of  fish  on  the  propagation  and  scattering  of  sound  in  the  ocean  and  others  have  been  conducted  in 
connection  with  studies  on  fish  sound  production  and  hearing. 

The  first  experimental  work  was  conducted  by  Coate  [36].  In  a short  series  of  measurements  on  a dead 
20-cm  long  crappie,  he  determined  that  a spherical  air  bubble  in  water  was  a reasonable  model  for  resonant 
scattering  from  a swimbladder,  provided  that  an  ariaitrary  loss  temi  was  included.  The  magnitude  and  form  of 
this  loss  term  were  left  unspecified. 

Batzler  and  Pickwell  [15]  made  measurements  on  several  live  goldfish  and  anchovies,  which  were  6-  to  11- 
cm  long.  The  experimental  resonant  frequencies  were  found  to  agree  well  with  those  predicted  by  equation  1- 
14  if  a value  of  5 x 10^  dynes/cm^  were  used  for  p,,.  However,  measured  Q's  were  in  the  range  of  3 to  5, 
which  are  lower  than  predicted  by  Andreeva. 
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McCartney  arxJ  Stubbs  [37]  made  measurements  on  live  cod,  saithe,  and  pollack,  wnich  were  30-  to  55-crri 
long.  The  experimental  resonant  frequencies  were  found  to  be  40  to  100  percent  higher  than  those  predicted 
for  a free  bubble  and  about  1 5 to  60  percent  higher  than  those  predicted  by  Andreeva.  The  differences  were 
attributed  to  the  ruxispherical  shapes  of  the  swimbladders  and  to  tension  in  swimbladder  tissues  Again 
measured  Q's  were  lower  than  those  predicted  by  Andreeva,  having  values  of  about  1 to  3,5. 

Sarxl  and  Hawkins  [33],  using  McCartney  and  Stubb  s measurement  system,  made  measurements  on  live 
cod  which  were  13-  to  50-cm  long.  When  the  fish  were  at  a depth  at  which  they  had  been  allowed  to  adapt  for 
several  days,  the  experimental  resonant  frequencies  were  found  to  be  from  about  25  to  600  percent  higher 
than  those  predicted  for  a free  bubble  and  the  measured  Q's  were  approximately  1.  The  ratio  of  the 
expenmental  resonant  frequency  to  that  of  a free  bubble  was  found  to  generally  increase  with  increasing  fish 
size.  Sand  and  Hawkins  suggest  that  the  cod  maintained  a tension  in  the  swimbladder  wall,  which  raised  the 
experimental  resonant  frequencies. 

Sand  and  Hawkins  compared  their  experimental  resonant  frequencies  to  those  McCartney  and  Stubbs 
obtained  from  cod  and  found  that  both  sets  agreed.  Sand  and  Hawkins  state  that  McCartney  and  Stubbs  did 
not  obtain  comparable  ratios  of  resonant  frequencies  because  they  misinterpreted  their  results 

Sundnes  and  Sand  [38],  also  using  McCartney  and  Stubbs’  measurement  system,  made  measurements  on 
live  chart,  which  were  probably  about  40-cm  long.  The  experimental  resonant  frequencies  were  found  to  be 
about  40  piercent  higher  than  predicted  for  a free  bubble.  About  three-fourths  of  the  difference  was  attributed  to 
the  shape  of  the  swimbladders  and  the  remainder  attributed  to  tissue  effects. 

Tavolga  [39]  and  Popper  [40]  have  examined  swimbladder  resonance  in  connection  with  fish  physiology 
studies.  In  both  cases,  probes  were  placed  inside  the  fish.  Tavolga  placed  a microphone  aganis*  the 
swimbladder  and  Popper  placed  a pressure  probe  inside  the  swimbladder,  puncturing  the  swimbladder  v^u.i  m 
the  process.  Tavolga  obtained  experimental  resonant  frequencies  which  were  two  to  three  times  higher  than 
predicted  for  a free  bubble.  However,  the  same  results  were  obtained  for  balloons,  whereas  others  have  found 
that  balloons  agree  quite  well  with  free  bubbles  [15, 37].  Tavolga  does  not  give  any  values  for  Q.  but  states  that 
it  is  quite  low.  Popper  found  the  pressure  response  inside  the  swimbladder  to  be  essentially  flat  and  therefore 
did  not  obtain  any  resonant  frequencies.  It  is  quite  possible  that  Tavolga's  and  Popper  s results  are  different 
from  those  discussed  above  because  of  the  influence  cf  the  probes  placed  inside  the  bodies  of  the  fish 

Love  [41]  conducted  an  experimental  study  which  was  related  to  swimbladder  resonance.  In  this  study,  a 
comparison  was  made  between  scattering  strengths  measured  acoustically  and  scattering  strengths 
calculated  from  coincidentally  collected  fish  distribution  data.  The  calculations  utilized  Andreeva's  equations  It 
was  noted  that  closer  agreement  between  calculated  and  measured  scattering  strengths  could  be  obtained 
near  resonance  if  lower  Q's  were  utilized  Although  this  study  does  not  present  any  direct  evidence  about  the 
proper  values  for  cot,  or  Q,  it  is  mentioned  here  because  its  implications  led  to  the  investigation  which  is  the 
subject  of  this  report. 

The  available  experimental  data  indicate  that  in  most  cases  the  resonant  frequency  of  a small  air  bubble 
gives  a reasonable  estimate  of  the  resonant  frequency  of  a fish  swimbladder,  especially  if  the  frequency  is 
adjusted  to  account  for  swimbladder  shape.  However,  the  resonant  frequencies  of  large  cod  were  found  to  be 
several  times  higher  than  those  predicted  for  air  bubbles.  In  addition,  in  ail  cases,  the  experimental  Q's  were 
lower  than  those  predicted  for  a swimbladder  fish  by  Andreeva 
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CHAPTER  II 


A NEW  MODEL 


The  Model 

The  purpose  of  this  report  is  to  Improve  the  equations  which  are  presently  utilized  to  predict  the 
resonant  frequerKy  and  acoustic  cross  section  of  a swimbladder  fish,  in  order  to  facilitate  the  correlation  of 
acoustic  and  biological  volume  reverberation  data.  A new  model  for  a swimbladder  fish  is  proposed  arxj  the 
appropriate  equations  developed  for  it.  The  model  consists  of  a small  spherical  shell,  enclosing  an  air  cavity, 
in  water.  The  shell  is  chosen  to  be  a viscous,  heat-conducting  Newtonian  fluid,  with  the  physical  properties  of 
fish  flesh,  and  the  interface  between  the  shell  and  the  cavity  supports  a surface  tension. 

The  shell  is  insonified  by  a harmonic  plane  compressional  wave  whose  wavelength  is  large  compared  to 
the  shell  diameter.  For  convenience,  the  center  of  the  shell  will  be  taken  as  the  origin  of  the  coordinate  system 
and  the  wave  will  travel  in  the  positive  direction  along  the  z-axis.  Spherical  coordinates  (r,  9,  rp)  will  be  used 
to  describe  the  field,  where  r is  the  distance  from  the  origin.  6 the  polar  angle  and  rp  the  azimuthal  angle. 
Since  the  problem  is  axisymmetric.  8/S(p  = 0. 

Acoustically,  fish  flesh  may  be  considered  to  resemble  soft  rubber.  In  such  a substance,  shear  waves  are 
much  less  important  than  transverse  waves  and  the  substance  can  be  closely  approximated  by  a viscous  fluid. 
Modelling  fish  flesh  as  a viscous  liquid  has  advantages  over  modelling  it  as  an  elastic  solid  with  a complex  shear 
modulus.  Specifically,  exact  equations  (Navier-Stokes)  exist  for  the  motion  of  viscous  fluids,  whereas  a 
phenomenological  approach  is  required  if  a complex  shear  modulus  is  utilized. 

The  body  of  a fish  which  surrounds  its  swimbladder  will  cause  an  increase  in  stiffness  over  that  of  a free 
bubble.  It  is  assumed  that  most  of  that  increased  stiffness  in  concentrated  in  the  swimbladder  wall  and  may  be 
modelled  by  a surface  tension  at  that  interface.  Unlike  an  elastic  modulus,  swimbladder  tension  may  be  under 
the  fish’s  control  and  could  account  for  otherwise  unexpected  variations  in  resonant  frequency. 


Limits  and  Physical  Properties 

Several  limits  will  be  placed  on  this  model.  These  include  limits  on  applicable  depth,  frequency,  and  size 
ranges.  The  depth  range  of  interest  is  from  the  surface  to  1,000  m,  which  corresponds  to  an  ambient 
pressure  of  10®  to  10®  dynes/criF.  The  frequency  range  is  100  Hz  to  40  kHz,  which  corresponds  to  a circular 
frequency,  co,  of  2n  x ICF  to  8n  x 10^  rad/sec.  The  fish  size  range  is  1 cm  to  1 m.  This  roughly  corresponds 
to  an  inner  shell  radius,  a,  of  1 0 ' to  5 cm.  (Appendix  A contains  a discussion  of  swimbladder  volumes.)  The 
ratio  of  outer  shell  radius,  b,  to  a is  2.5  S b/a  S 6 for  small  fish  and  2.5  < b/a  < 3.2  for  large  fish.  Two  further 
limitations  will  be  that  6 x 10^  cm/sec  < wa  < 2.4  x lO''  cm/sec  and  that  a)b  < 7 x 10^  cm/sec. 

In  addition  to  swimbladder  size,  the  surface  tension  at  the  swimbladder  wall  must  be  specified  For  small 
fish, 

1CF  dyne/cm  < s < 10®  dyne/cm 

and  for  large  fish, 

10^  dyne/cm  < s < 10®  dyne/cm. 

The  rationale  for  this  choice  of  ranges  is  discussed  in  Appendix  A. 

Several  other  parameters  must  also  be  specified.  These  include  the  velocity  of  a compressional  wave 
(sound  velocity),  c;  density.  Po;  specific  heat  at  constant  pressure,  Cp;  ratio  of  specific  heats,  y:  thermal 
conductivity,  k;  shear  viscosity,  ri,:  and  bulk  viscosity  rip;  for  air,  sea  water,  and  fish  flesh.  Subscripts  a.  w, 
and  f will  be  utilized  to  indicate  the  properties  of  air,  sea  water,  and  fish  flesh,  respectively.  For  convenience, 
a viscosity  parameter,  is  defined  as 

^ = 30,. + nw  (11-1) 

The  parameters  are  listed  in  Table  1,  for  a temperature  of  10°C  and,  unless  specified,  a pressure  of  one 
atmosphere  (10®  dyne/cm^).  The  properties  of  air  were  obtained  from  references  42  and  43  and  those  ot 
sea  water  from  reference  44.  The  properties  of  fish  flesh  are  discussed  in  Appendix  A 
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Table  1 — Physical  Properties 


Air 

Sea  Water 

Fish  Flesh 

c,  cm/sec 

3.3  X 

10“ 

1.50  X lO'’ 

1.55  X 10^ 

Po,  gm/cm^ 

1.3  X 
1.3  X 

10  ^ (at  1 atm) 

10  ’ (at  100  atm) 

1.026 

1.050 

Y 

1.40 

1.01 

1.01 

cal 

gm  °C 

0.24 

0.93 

0.89 

K, 

5.5  X 

10  5 

1.34  X 10  3 

1.32  X 10  3 

cm  sec  u 

ris,  poise 
poise 

1.8  X 

10  “ 

1.4  X 10  2 

1 to  10“ 

Preliminary  analysis  of  the  new  model  indicated  that  considerable  simplification  resulted,  with  virtually 
no  loss  of  accuracy,  if  several  approximations  were  made  concerning  the  physical  properties.  The  first 
approximations  were  that  y*  = Yt  = 1 . Secondly,  since  k,>ok3,  it  was  assumed  that  any  heat  generated 
in  the  air  is  rapidly  conducted  away,  so  that  the  temperature  in  fish  flesh  and  water  is  constant.  Thirdly, 
since  ris,»ris^»r|sg,  it  was  assumed  that  r|s„  = = 0. 

One  further  assumption  which  greatly  simplifies  the  problem  is  made  in  this  model.  As  discussed  in 
Chapter  I,  only  the  fundamental,  or  zeroth  mode  of  oscillation  will  be  considered. 


Formulation  of  Equations 

The  first  step  in  the  determination  of  the  resonant  frequency  and  acoustic  cross  section  of  the  new 
model  is  the  determination  of  the  oroper  wave  equations  in  water,  fish  flesh,  and  air.  The  wave  equations 
are  determined  from  the  basic  equations  of  motion  as  given  by  Hunt  [45],  generally  following  the  method 
discussed  by  Epstein  and  Carhart  [46].  The  basic  equations  are  the  continuity  equation, 

|£  + V • (pu)  = 0 ; (11-2) 


the  equation  of  conservation  of  momentum, 

p = - p(u  • v)u  - VP  -I-  ^v(V  u)  - ; ("'3) 

o t 

the  equation  of  conservation  of  energy, 

pc,  ( + u • VT ) + V u+V  q-<p,  = 0;  (11-4) 

and  the  equation  of  state, 

p = P(P,T)  : (11-5) 

where  p is  density,  t is  time,  u is  a velocity  vector,  P is  pressure,  c,  is  specific  heat  at  constant  volume,  T is 
absolute  temperature,  (3  is  the  coefficient  of  thermal  expansion,  "q  is  a heat  conduction  flux  vector, 

q = -kVT  , ("-6) 

and  (p,  is  a viscous  dissipation  function. 


and 


Equations  11-2  through  11-5  are  linearized  by  assuming  that 
u = 0 + u,  . 
p = Po  + Pi  . 

T = To  + T,  . 


(11-7) 
(11-8) 
(11-9) 

P(P.T)  = Po(Po.To)  + Pi  • 
where  the  subscripts  0 and  1 refer  to  the  average  and  perturbation  values  of  the  parameters,  respectively 
The  first-order  equations  are 


^ + Po  '^^l)  ~ ® • 


3 t 


Po-^  = - 7P,  + EV(V  u,)  - ii.7x(Vxu,)  , 


(11-11) 

(11-12) 


PoCv  ^ + Po-%  — V ui  - = 0 . 


(11-13) 


since  cp,  is  of  second  order  [45],  and 

^ (it)/' 


(11-14) 


It  has  been  assumed  that  T„  = T,,  = 0 and  that  Yi  = Y.  = 1 Hence,  in  fish  flesh,  equations  ll-l  1 through 
11-14  are 


dp, 


^ + Po,(7  • u,)  = 0 , 

Po,-|^  = - VP,  + EV(7  u,)  - n.,vx(7xa,)  , 


(11-15) 


and 


K _ 1 (-A£s—\ 


, ...  (II-16) 

" (■^)t 

These  equations  are  the  same  in  water,  with  the  exception  that  the  viscosity  is  set  to  zero.  The  isothermal 
compressibility,  Kt,  is  defined  as 

(11-18) 

(11-19) 
(11-20) 

(i:-21) 


and  the  adiabatic  sound  velocity,  c,  as 

C2 


PoKt 

Since  Y,  = Yw  = T equation  11-17  is 


Pi,  = 


_ Pi 


0,2 


and  equation  11-15  is 


^ + Po,c,2(V-u,)  = 0 . 

According  to  Helmholtz's  theorem  [47],  the  u,  vector  field  can  be  lopresented  in  terms  of  a vector 
potential  A and  a scalar  potential  Q such  that 


1 


‘ 

1 

u,  = V X A - VO  . 

1 

(II-22) 

with  V A = 0 . 

(II-23) 

i 

so  that  V u,  = - V^O  , 

(II-24) 

Ik 

i 

and  VxD,  = - WA  . 

(II-25) 

Thus,  equations  11-21  and  11-16  become 

t 

r 

Po,c,^V^n  = -|^ 

(II-26) 

f 

and 

b 

i 

Po,  -|j-{V  X A)  - ns,V  X (V^A)  = - VP,  - ^V(V^Q)  + Po,  • 

(II-27) 

f 

f 

1 

If  the  curl  of  equation  11-27  is  taken, 

V X VP,  = V X VO  = V X V(V20)  = 0 . 

(II-28) 

f 

j 

Thus, 

^Is 

Po,  -|y(V  X A)  - ns,V  X {V^A)  = 0 , 

(II-29) 

* 

or 

ns,v^A  - p„,  = 0 . 

(I1-30) 

Therefore, 

VP,  + ^V(V20)  - Po,-|j-(VO)  = 0 . 

(I'-SI) 

1 

Then  substituting  equation  11-26  into  11-31  and  taking  the  divergence  yields: 

V’r  i ^ v=-  'i  - — ^ 

^ ^ Po,c,^^  V at  ; C,^  at= 

(11-32) 

/l' > 

Equations  11-30  and  11-32  are  the  wave  equations  in  the  fish  flesh.  In  water,  where  the  viscosity  is  zero,  the 

wave  equation  is 

f' 

: % 
ii 

n.p  _ 9^'"'  - 0 

^ c„2  aT’ 

(11-33) 

i' 

The  harmonic  time  dependence  is  now  introduced,  such  that 

p,  = Pje  . 

(11-34) 

T,  = T2  e , 

(11-35) 

’; 

and 

/ 

>1 

II 

>1 

(D 

1 

(11-36) 

Then  equations  11-30, 11-32,  and  11-33  are,  after  eliminating  the  numerical  subscripts. 

< 

1 

,> 

10 

j 

n„V2A  + iu)po,A  = 0 . 


V Po,Ct'  / 


V^P  + -~P  = 0 . 
C,^ 


y^p  + -iii-p  = 0 . 


The  coefficient  of  thermal  expansion.  3,  is  defined  as 


B = _ _L 

^ On  ' dT  /P 


^ Po  '^T  /P' 

so  that,  in  air,  equation  11-14  is 


P.3  = 


— — D B T 

^TT  PoaPa*.  • 


In  addition,  equations  11-1 1 through  11-13  for  air  are,  after  substitution  for  u: 
- Po^V^O  = 0 . 

Po3-|j(V  X A)  = Po,  '^(VO)  - VP,  , 


9 T,  (Cp,  c„,) 


Poa^va  0 I Po, 


V=Q  - KpV^T,  = 0 . 


Taking  the  curl  of  equation  11-43  and  utilizing  equation  11-28  indicates  that 


so  that 


V X A = 0 . 


VP,  - Po,  -^(VO)  = 0 . 


Then,  substitution  of  equation  11-41  into  11-42  yields: 


Hence,  taking  the  divergence  of  equation  11-46  and  substitution  of  equation  11-47  yields; 

In  addition,  substitution  of  equation  11-47  into  11-44  yields; 

Po.c„.  Ill  - - K,V^T,  = 0 . 


►'Oa-Pa  0t  P,C,2  at  “ ' 

Introduction  of  the  harmonic  time  dependence  into  equations  11-48  and  11-49  yields,  after  eliminating  the 
numerical  subscripts; 


V«p  + i^^P  - (*)2po,PaT  = 0 


(11-50) 


--  - • A..  • ■ ^ 


11 


'ft  i 
11  ! 


Since  k,,  and  k2,  are  never  equal,  the  general  solution  of  equation  11-53  is 
P = Hj,  + 412  . 

where  41,  and  4)2  satisfy  the  equations 

(V2  k,,2)  4<i  = 0 

and 

(V*  + k2,2)  4*2  = 0 . 

Equations  11-37  through  11-39  can  be  written  in  a form  similar  to  equations  11-62  and  11-63.  Thus 


(11-62) 

(11-63) 


(72  + k3,2)  A = 0 , 

(11-64) 

(V2  -f  k2,2)  P = 0 , 

(11-65) 

and 

(V2  + k2,2)  p = 0 . 

(11-66) 

where 

k 2 _ '“Poi 
''31  n.  ' 

(11-67) 

c,2  Po.c,2  ) ' 


k 2-J£. 

""  "c^ 


Comparisons  of  the  equations  for  air  and  water  indicate  that  4^2  represents  a compressional  wave.  Hence,  4) , 
represents  a thermal  wave. 

Equation  11-64  can  be  transformed  into  a scalar  equation.  It  can  be  shown  that  a vector  f exists  such  that 
[48] 

u = V^r  (11-70) 


u 

% » 

■ [■ 
f I 


s = V X r . 

Since  the  problem  is  axisymmetric,  8/3  cp  = u,  = 0,  which  implies  that  r,  = 0.  Thus 

Vxr.r(0)  + e(0)  + 5;-;-[-^(-r.,-.^]  . 

so  that 

A = . 

where?,  5,  and  (p  are  unit  vectors.  Therefore 
7^.$  [VA,-^^]. 


(11-71) 

(11-72) 

(11-73) 
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and  equation  11-64  can  be  written  as 


The  plane  wave  impinging  on  the  shell  will  be  represented  by 

P»,  = Ae*2w^  . (11-76) 

where  P„,  is  the  perturbation  pressure  of  the  wave  in  water,  and  A is  the  pressure  amplitude  P,  can  be 
expanded  in  spherical  waves  as 

at 

p«,  = A J (2  e + l)i^jy(kjj)  Pg  (cos  0)  . 

K-O 


where  (!  is  the  mode  number,  )(,(k2^r)  is  a spherical  Bessel  function  and  P^  (cos  0)  is  a Legendre  polynomial 
The  incidence  of  the  wave  upon  the  shell  gives  rise  to  five  additional  waves,  which  are  represented  by  the 
solutions  of  equations  11-62,  11-63,  11-65,  11-66,  and  11-75.  These  equations  can  be  solved  by  standard 
separation  of  variables  techniques.  The  solutions  are: 

4>1  = ^ Daejg(kiaO  Pe(cos0)  , 

(11-78) 

X 

4*2  = X B,gig(k2/)  Pk(cos  0)  , 

(?-0 

(11-79) 

X 

Pi  = ^ (kj.r)  -1-  E,  gng(k2,r)]  Pg(cos0)  , 

(11-80) 

X 

P8(cos  0)  , 

(11-81) 

and 

"*■  G,gng(k3,r)]  P'£(cos  0)  , 

(11-82) 

where  B,  D,  E,  F,  and  G represent  the  amplitudes  of  the  waves,  ri  Jk2,r)  is  a spherical  Neurriann  function, 
and  h ^ (kj^r)  is  a spherical  Hankel  function  of  the  first  kind.  P*,^  indicates  the  scattered  compressional  wave  in 
water,  so  that 


Pw 


-t^  P. 


(11-83) 


One  of  the  assumptions  of  this  model  is  that  only  the  fundamental  mode  contributes  to  the  scattering.  Thus 
only  the  6 =0  mode  is  considered.  Examination  of  equation  11-82  shows  that  for  C = 0, 


A„  = 0 . 

The  remaining  equations  are  now  written  for  C = 0: 

(II  84) 

Pw  = Ajo(k2„r)  + B*ho(k2„r)  , 

(11-85) 

P|  ~ B|jo(k2|r)  -I-  EfnQ(k2|r)  , 

(11-86) 

ij)!  = DJo(k,/)  , 

(11-87) 
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and 


>^2  = Bjo(kj.r)  , 


(11-88) 


I 


r 

« 


» 


where  the  numerical  subscripts  have  been  eliminated  from  the  coefficients  for  convenience. 
The  scattenng  cross  section  of  the  shell  is 

o = ®,/l,  . 

where  9,  is  the  average  scattered  power  and  I, the  incident  intensity. 

where  v is  the  area  of  an  element  of  a sphere  (49). 

V = Rv  , 
dv  = r^dx  . 

where  n is  the  normal  to  the  area  and  x is  the  solid  angle.  Also, 

I,  = 4 “«■ 

so  that  it  is  necessary  to  determine  u,,^  and  u„  in  terms  of  P,,  and  P*,,  respectively  [49]. 

In  the  fundamental  mode  u = - Vfi,  so  that  utilizing  equation  11-31  for  water, 


(11-89) 

(11-90) 

(11-91) 

(11-92) 

(11-93) 

(11-94) 


dependence  yields; 


9 r 

ap,, 

0Z 

Thus,  from  equations  11-76  and  11-97, 
2Po»Cw 

From  equation  11-81, 

P«,  = B,.ho(k2«r), 


Introducing  the  harmonic  time 


so  that 


and 


= — P, 
* wPo* 


wPo» 


~ “Pow 


where 


e*2w' 


(11-95) 

(11-96) 


(11-97) 


(11-98) 

(11-99) 

(11-100) 
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J 


1 


1 


so  that  for  r — oc , 


ho'(k2^r)  = 


(11-101) 


Thus,  as  r — oc, 


P * M = 

"»  wpo.k2. 


(11-102) 


2nB,B,‘ 


(11-103) 


Therefore, 


4nlB»|2 


(11-104) 


and  in  order  to  determine  o,  it  is  necessary  to  solve  equations  11-85  through  11-88  for  (B„/A), 


Boundary  Conditions 

Equations  11-85  through  11-88  contain  five  unknown  coefficients.  Thus,  five  boundary  conditions  are 
required  for  the  solution  of  this  set  of  equations.  These  conditions  are  the  continuity  of  normal  velocity,  u„  and 
normal  stress,  t„,  at  both  interfaces  and  continuity  of  temperature  at  the  inner  interlace.  The  first  step  in 
solving  for  B«/A,  therefore,  is  to  obtain  u„  and  T in  terms  of  P*,  P,,  if),  and  ifij. 

Equation  11-50  gives  T^  in  terms  of  Pg.  Substitution  of  equations  11-58  and  11-60  through  11-63  into  equation 
11-50  leads  to 


E7J7  I [v. 


(11-105) 


Urw 's  given  in  terms  of  P„  by  equation  11-96.  Following  the  same  procedure  to  obtain  u„,  equation  11-31  may 
be  written  as 


p,  + - po.  ^ = 


(11-106) 


Then  substitution  of  equation  11-26  into  11-106  yields,  with  intrcxJuction  of  the  harmonic  time  dependence: 

^ r 1 _ 1 D ... 


(11-107) 


Similarly,  from  equation  11-46  and  11-61 : 


(11-108) 


i_  r,  _ 1 ^ 

topo,  I PofCi^J  ar 


(11-109) 


= - _L  r 

'*  WP03  I 9r  9r  j 


(11-110) 
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The  normal  stress  in  fish  flesh  is  [45] 


= - P,  + . 


(11-111) 


so  that  utilization  of  equation  11-109  yields: 


T = -P  -Ji  fi 

(DPo,  L Po,c,U  dr^ 


(11-112) 


T = - P 


(11-113) 


— ~ Pa 


(11-114) 


At  the  air  fish-flesh  interface  the  surface  tension  must  be  included  in  the  stress  equation.  For  a bubble  in 
water  [50], 

P«,  = Poo,+-^.  (11-115) 


(II-115) 


where  P and  Pout  are  the  total  pressures  inside  and  outside  the  bubble  and  R is  the  radius  of  the  bubble  at  any 
instant.  If  g represents  the  small  changes  in  the  radius  of  the  bubble  and  a is  the  average  radius,  such  that 


R = a + g , 


(11-116) 


where  g<<a,  then 


R a(1  g/a)  a 


- — (1  - g/a)  , 


(II-117) 


so  that 


p,„  = Pou.  (1  - g/a)  . 


(11-118) 


Since  P,nand  Pg^,  are  total  pressures,  equation  11-1 18  can  be  linearized  utilizing  equation  11-10.  Thus,  the  first 
order  equation  is 


P = p + 

' /Mill  ' *11  ' 


outi  ' mi  ’ q2  * 


(11-119) 


Due  to  the  harmonic  time  dependence. 


IT  - -'“8  ■ 


(II-120) 


p _ p = 

^ oull  ' mi  (jjgZ 


(11-121) 


The  boundary  conditions  are: 


BC  1 u,^  = u„  at  r = b , 


(II-122) 


BC  2 at  r = b , 


(11-123) 


'■  Tt**' 


BC  3 

Um  = 

at  r = a , 

(11-124) 

2iUraS 

BC  4 

^rrt  — 

-1^  ^-3. 

(11-125) 

BC  5 

Ta  = 0 

at  r=a. 

(11-126) 

Equations  11-122  through  II-126  have  been  labeled  BC  1 through  BC  5 tor  later  convenience.  Substitution  ot 
equations  11-105,  11-109,  11-110,  11-96,  11-112,  11-113,  11-114,  and  11-85  through  11-88  into  equations  11-122 
through  11-126  yields: 

BC  1,  (-^)  [ AkjJo'(k2«b)  + B.kj„ho'(k2.b)  ] 

= ^-p  ^ 1 ~ p ^2  ^1^  BikjJg  (kjjb)  E,k2,no  (kp|b)  j , (11-127) 

BC  2,  Ajo(k2^b)  + B„ho(k2„b)  = B,jo(l<o,b)  -h  E,no(k2,b) 

^ [ ®'W(k.b)  . E,k,^ng"(k2,b)  ] , (11-128) 

BC  3,  ( (l  - -^^2  ) [ B,k2,Jo'(k2,a)  + E,k2,no'(k2,a)  ] 

= (-^  )[  D,k,Jo'(k„a)  + B,k2,jo'{k2aa)  ] : (11-129) 

BC4,  B,jo(k2,a)  + E,ng(k2,a)  + (-iL  ) ( i - [ B,k2,^jo''(k2,a)  + E,k2,^no"(k2,a)  ] 

= Daio(k,aa)  + Bajo(k2aa)  ^ ) [ Dak,Jo'(k„a)  - Bak2,io’(k2,a)  ] ; (|M30) 
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CHAPTER  III 

SOLUTION 


The  first  step  in  the  solution  of  equations  11-127  through  11-131  for(B„  /A)  is  to  solve  11-131  for  D,  and 
substitute  the  result  into  equations  11-129  and  11-130.  For  conciseness  the  boundary  conditions  will  now  be 
written  in  a tabular  form: 


BC  B.  B, 

1 S„  S,j 

2 Sj,  S22 

3 0 832 

4 0 S42 

The  S,’s  are  taken  such  that 


S,iB^  -t-  S,2Bf  -t-  S,3E|  -t-  S,4B,  — q,A 


The  S,|’s  and  q’s  are  as  follows: 

521  = — ho(k2„b)  , 

522  = jo(l^2tb)  ( pp  ) io”(l^?ib)  ■ 

523  = no(k2,b)  + ( pg  ) rio"(k2,b)  , 

^2  = jo  (k2wt>)  . 

" ( Po,Sk2,)  ’ 

(Ya-1)  k,.jo'(k,.a)jo(k23a) 
(Va  jo(k„a) 

542  = jo(k2,a)  + ( p^q2  ) io"(k2,a)  . 

543  = no(k2,a)  + ( no''(k2,a)  . 


k2aio'(k2aa) 


A 

a, 

02 

0 

0 


(MI-1) 

(III-2) 

(MI-3) 

(III-4) 

(III-5) 

(III-6) 

(III-7) 

(III-8) 

(III-9) 

(111-10) 

(111-11) 

(111-12) 

(111-13) 

(111-14) 
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(111-15) 


S44  ~ 


(Ya-1)jo(k;,a) 

/ _iP03Ca^Cp,\ 


r . / 2s  \ / k,Jo'(k,.a)  \i 
L ^ V co^a^Po^  / V Jo(k'.a^  / . 


+ (l||^J  io'(k^.a)  . 


where  equation  11-68  has  been  used  for  kj,^. 

The  solution  of  the  boundary  conditions  for  is 

B„  a,U  + OjW 
“S"  S„U  + SjiW 

where 


and 


U — 244(322833  S23S32)  + 334(323842  22384 


W — 844(343832  212833)  + 834(812843  813842) 


(111-16) 


(111-17) 


(111-18) 


This  solution  is  easily  arrived  at  after  several  pages  of  substitutions. 

Equation  111-16  could  be  solved  numerically  with  a computer  by  choosing  values  for  the  various  physical 
parameters  and  calculating  the  3, , 8 an'^  a,'s.  However,  this  will  not  be  done  because  the  purpose  of  this 
report  is  to  obtain  simplified  equations  for  Up  and  o which  can  be  solved  without  resorting  to  use  of  a 
computer.  The  simplified  solution  is  based  on  the  ranges  of  the  physical  properties  and  the  limits  of  the 
variables  given  in  Chapter  II.  It  is  obtained  by  assuming  that  I ki^a  | islargeandthatk2„b,  k2,b,  andk23aare 
small  and  accepting  any  errors  which  are  less  than  ten  percent.  Order  of  magnitude  comparisons  are  then 
made  between  various  terms  and  any  term  which  is  always  less  than  ten  percent  of  another  term  is 
negl'.cted.  This  simplification  process  is  shown  in  Appendix  B.  The  result  of  this  process  is; 


= - (1  + A + iX) 


1 + 


[ 


+ I 


3Po,C«(Ya  - 1) 
Po*“=a^ 

CwPof 


(^)[(4^)('-Spj7^)-’-(W^)]!' 


(111-19) 


where 


and 


^ ^ / / 11Po,C,2  _ \ 

V Poi^Ci-as  / V 9po„c*2  ) 

/ 2a)^b3  \ r / pp,c,=  _ , \ _ / Poi^  ^ n \ 1 

^ - V 3p;rc,^i^ ) L > VpowC.^  4 Hb3j  J 


(111-20) 


(111-21) 


Equation  111-19  can  now  be  substituted  into  11-104  to  yield  o.  However,  before  doing  this,  equation  111-19 
will  be  manipulated  so  that  the  final  results  can  be  easily  compared  to  those  of  other  researchers.  Thus,  the 
resonant  frequency  is  defermined  by  setting  the  imaginary  part  of  the  denominator  of  equation  III-1 9 equal 


(111-22) 


( 3P°aC.==  \ /i  _ 2S  1 _/ 

V w^azpo, ) 1.  3po,Ca2a  ) I 


\ 

}po,Ci^a2  } 


and  solving  for  oj^a^  yields  the  resonant  frequency  of  the  new  model: 


3p0aC,2 


(l — 

V 3po3Ca2a 


3po,2c,2a2 


Now  [42], 

PoaCa^  = YaPo 

and  from  equation  11-115, 


P — P -t-  — p + 

^Oa  - r^Of  + g - r^o„  + g 


Therefore,  substitution  of  equations  III-24  and  III-25  into  III-23  yields: 


tOo^a^  = 


+ — (3Ya  - 1) 
Po(  Pot® 


3po,2C,2a2 


The  damping  factor,  H,  will  be  defined  by 


Po.<»a  , _2i 
(UH  Pof^w  Poia>< 


PofCw  Po,a)a2 

3(Ya  - 1)  / WKa 


V 2PoaCpa  } \ 


Po,co2a3 


Thus,  following  Devin  [19], 


JL  = _L  -L  -L  = H 

rtrao  '■vij  n,(, 


where 


w _ t*1oPoi<'w 
r'rad  “ ^ ' 


a)2po«a 


H.,.  . , 


U _ '"O'* 

3(Ya  - 1) 


ioa  /2po,Cp  ^ 2s  \ ' 

_ 1)  V COKa  / V Po,t«)^a3  / 


(III-23) 


(1I1-24) 


(III-25) 


(III-26) 


(III-27) 


(III-28) 


(III-29) 


(III-30) 


(1II-31) 


As  before,  at  co  = Wq,  H = Q.  Thus,  substitution  of  equations  111-26  and  111-27  into  111-19  yields: 


CHAPTER  IV 

DISCUSSION 


Results 

The  equations  developed  for  (Jq,  o,  and  H In  Chapter  III  for  the  new  model  can  be  considered  to  be 
accurate  within  ten  percent  only  for  the  parameter  ranges  given  in  Chapter  II  The  ranges  of  lor  which 
exact  values  are  not  known,  and  s,  which  may  be  varied,  were  extended  beyond  the  limits  of  present 
available  data  However,  the  largest  values  of  a,  w.  and  Po„  used  in  this  model  are  not  beyond  the  range  of 
possible  values.  These  limitations  weie  imposed  by  the  assumption  that  kj^a  is  small  This  required  that,  for 
the  error  caused  by  using  just  the  first  term  of  the  spherical  Bessel  function  expansions  to  be  less  than  ten 
percent, 

(joa  < 2,4  X ICt*  cm/sec. 

The  limitation  of  k^^^a  being  small  is  a factor  in  all  the  previous  models,  but  in  most  cases  it  is  not  mentioned 
or,  if  it  is  mention^,  no  limiting  value  is  given. 

The  range  of  ^ utilized  in  the  solution  of  the  new  model  was 

1 poise  < ^ <10“  poise. 

However,  a more  likely  range  of  as  shown  in  Appendix  A,  is 
50  poise  <^<2x10^  poise. 

This  narrower  range  will  be  utilized  in  the  remainder  of  the  discussion  because  it  results  in  a simplification  of 
equations  111-26  and  111-34  without,  in  all  probability,  affecting  the  results.  Equations  111-26  and  111-34  simplify 


/ Pow  \ 
VPo,  / 


o 

3 

1 (I)3HJ 

V 0)3  ' j J 

respectively,  when  the  upper  limit  of  ^ is  reduced  to  2 x 10^  poise.  A and  X are  negligible  compared  to  1 tor 
2 X W poise.  A and  X were  the  only  remaining  terms  which  contained  b.  Thus,  for  the  a/b  range  given 
in  Chapter  II,  b is  unimportant  for  ^ < 2 x 10^  poise. 

Figure  2 gives  cooa  as  a function  of  depth  and  surface  tension,  calculated  from  equation  IV-1.  s has 
essentially  no  effect  on  for  s/a  < 10^  dyne/cnr’  at  any  depth  and  for  s/a  < 10^  dyne/cm^  for  depths 
below  100  m.  For  s/a  values  around  5 x io^  dyne/crrF,  pressure  has  essentially  no  effect  on  (Dpa  for  depths 
less  than  100  m and  at  a depth  of  1000  m,  cOoa  is  increased  by  only  50  percent  over  its  value  near  the 
surface 

At  resonance, 

/ n \ 2 


o = 4na^O^ 


/_Pow_V 

V Po.  / 


Q<ad  ~ 


_ Pot^w 


Po„o)oa 


Q...  - . 

and 

^ uio'^a  / 2po,Cp.\  ^ f 

~ 3(y,  - 1)  V Ka  / \ 

Thus,  o is  proportional  to  O’  at  resonance 


2s 

Po,w<,*a3, 


Figure  3 shows  and  Q„s  as  functions  of  Woa.  increases  with  increasing  Wga  and  decreases  with 
increasing  ^/a.  Values  of 

2 X 10^  poise/cm  < ^/a  < 10^  poise/cm 

give  values  of  which  are  comparable  to  values  obtained  for  Q in  near-surface  swimbladder  resonance 
experiments  [15.  33,  36-38].  Q,^  decreases  with  increasing  Woa  and,  depending  upon  the  value  of  E/a,  may 
be  significant  for  all  Uga  considered. 

Figure  4 shows  Q^,  as  a function  of  CDga,  a,  s/a,  and  depth  D.  Any  individual  numbered  curve  indicates 
that  Q,^  increases  gradually  with  increasing  s/a.  The  set  of  numbered  curves  indicates  that  Q,,,  increases 
with  increasing  a and  D.  The  lettered  curves  are  for  the  case  where  the  effect  of  surface  tension  is 
insignificant,  and  also  indicate  that  Q*,  increases  with  increasing  a and  D.  A comparison  of  figures  3 and  4 
shows  that  Q#,  is  significant  only  at  depths  above  100  m and  then  only  it  ^/a  < 1CF  poise/cm. 

At  off-resonance  frequencies 


(IV-7) 


(IV-8) 

(IV-9) 


Thus,  at  frequencies  above  resonance,  may  be  the  dominant  damping  term,  even  for  relatively  high 
values  of  ^a.  Conversely,  at  frequencies  below  resonance,  Kao  may  be  the  least  important  damping  term 
The  variation  of  with  (cjb/w)  is  complicated  by  the  presence  of  the  surface  tension  term.  It  can  be  shown, 
utilizing  equation  IV-1,  that 


Po.Wo^aa 


S 1.31 


(IV-10) 


so  that,  for  the  purpose  of  establishing  a trend, 

H -H  2s  V’n  . (IV-11) 

~ \loJ  \ P^;u7a3/ 

Hence,  at  frequencies  above  resonance  and  at  frequencies  below  resonance  for  which  po,Ci)-a  »(2s/a), 
~ Q.h  (1V.12) 


H.  (^)  Q„ 

At  frequencies  below  resonance  for  which  (2s/a)>>  Po,(D^a2, 

H„  , (3)’(P^)  O.. 

Thus,  becomes  relatively  less  important  than  as  frequency  increases  and  conversely. 


(IV-13) 


Comparison  to  Free  Bubble 

The  validity  of  the  results  developed  for  the  new  model  must  be  determined.  This  will  be  done  by  first 
checking  that  the  equations  developed  for  the  new  model  approach  the  equations  for  a free  bubble  in  water 
as  a limiting  case  and  then  comparing  values  calculated  using  the  new  model  with  experimental  values. 

The  equations  obtained  in  Chapter  III  for  the  resonant  frequency  and  scattering  cross  section  of  the  new 
model  can  be  readily  compared  to  the  equations  given  for  a free  bubble.  Equations  1-6  and  1-7  give  tOb  and  o 
for  an  ideal  bubble,  neglecting  surface  tension.  Equation  111-26  reduces  exactly  to  equation  1-6  when  ^ = s = 
0.  If  thermal  losses  are  also  neglected,  equation  111-33  reduces  to  equation  1-7  when  ^ = s = 0,  except  for 
factors  of  (Pov,/Pot)  If  surface  tension  is  included,  a)o  for  an  ideal  bubble  is  given  by  equation  1-12  with 
ris,,  = 0.  Equation  111-26  reduces  exactly  to  equation  1-12  when  ^ = ri,^  = 0.  Thus,  when  viscous  and  thermal 
dissipation  effects  are  eliminated  from  the  new  model,  the  results  are  equivalent  to  results  obtained  for  an 
ideal  bubble. 

When  the  viscosity  of  water  is  considered  for  an  air  bubble  in  water,  u)p  is  given  by  equation  1-12.  A 
comparison  of  equation  111-26  with  1-12  shows  that  the  form  of  the  viscosity  factor  for  the  new  model  differs 
from  that  given  for  a bubble.  However,  as  was  noted  earlier,  equation  1-12  is  only  valid  for  small  viscosities, 
so  that  this  difference  is  not  surprising  (25).  Both  equations  1-12  and  111-26  indicate  that  decreases  as  the 
viscosity  increases,  which  should  be  the  case  for  a viscously  damped  system  [51] 
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FIGURE  4 THERMAL  DAMPING  FACTORS  FOR  THE  NEW  MODEL  OF  A SWIMBLADDER 
BEARING  FISH 
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The  damping  in  the  new  model  can  also  be  compared  to  that  for  a real  bubble  in  water  as  determined  by 
Devin  [19).  For  ease  in  comparison  with  Devin’s  results,  equations  IV-4,  IV-5,  and  l\/-6  will  be  rewritten  as: 


and 


Po><»>oa 

2^ 

(»)oPo,a2 


(IV- 14) 
(IV-15) 


6,h 


3<Ya  - 1)/  K, 

Wo'’a  V 2po.Cp.  ) V 


Po,u)o=a3  ; 


(IV- 16) 


Equation  IV-14  is  equivalent  to  equation  1-9,  except  for  the  factor  (Po,/Po,)  Devin  has  considered  only  the 
shear  viscosity,  so  that  equation  11-1  appears  to  indicate  that  equations  IV-15  and  I-10  differ  by  a factor  of 
Vi.  However,  and  r),  are  related  to  the  dilatafional  viscosity  Pa,  by  [45] 

Hb  = Hd  + ■ (IV-17) 

Therefore,  if  ria  is  neglected,  equation  ll-l  indicates  that  ^ = 2ri5,  Hence,  equations  IV-15  and  1-10  are 
equivalent.  Equations  IV- 16  and  1-1 1 are  also  equivalent,  except  for  the  surface  tension  factor  in  equation 
IV- 16.  Devin  considered  surface  tension,  but  since  the  surface  tension  of  an  air  bubble  in  water  is  only  74 
dynes/cm,  its  effect  is  insignificant  for  the  frequencies  and  bubble  sizes  of  present  interest.  Hence,  the  form 
of  the  damping  factor  at  resonance  for  the  new  model  is  essentially  equivalent  to  that  of  a real  bubble  in 
water. 


Comparison  to  Experimental  Data 

The  previous  comparisons  compared  the  spherical  model  to  spherical  bubbles.  In  order  to  compare  the 
model  to  experimental  measurements  made  on  swimbladder-bearing  fish,  it  is  obvious  from  figure  1 that  the 
effect  of  swimbladder  shape  should  be  included.  This  is  done  by  letting 

Woc  = ^“>0,  . (IV- 18) 

where  cut,  is  the  theoretical  resonant  frequency  for  the  spherical  model,  ? is  given  by  equation  1-25,  arxJ 
is  the  value  which  is  utilized  in  the  comparisons  to  the  experimental  data 

There  are  five  sets  of  experimental  data  to  which  calculations  based  on  the  new  model  can  be 
compared.  These  are  contained  in  references  15,  33,  36,  37,  and  38.  In  each  of  these,  the  measured 
resonant  frequency,  (Oq^,  and  the  calculated  Q are  either  given  directly  or  can  be  obtained  from  a curve  The 
data  in  reference  36,  which  was  collected  a decade  before  any  of  the  other  data,  have  variations  which 
require  some  interpretation  to  obtain  ut)„,  and  Q.  This  required  interpretation  makes  the  u)o„,  and  Q values 
obtained  from  reference  36  less  reliable  than  those  obtained  from  other  sources 

Comparisons  of  the  data  to  the  model  require  that  values  for  a and  e be  krrown.  In  some  cases,  a and  e 
were  measured  and  the  values  given.  In  the  other  cases,  it  was  possible  to  determine  a indirectly  from  other 
data  in  the  report.  When  values  of  e were  not  given,  average  values  were  obtained  from  other  sources 
Measurements  of  a and  e were  all  made  at  atmospheric  pressure  There  is  some  question  as  to  what  effect 
increasing  depth  has  on  swimbladder  volume  [e  g.  11.  12];  that  is,  whether  the  fish  retains  a constant 
swimbladder  volume  or  a constant  swimbladder  mass  (Boyle's  Law)  or  some  other,  intermediate  process 
(McCartney  and  Stubbs  [37]  assumed  an  intermediate  process,  due  to  the  effect  of  tension  in  the 
swimbladder  wall.)  It  is  probable  that  different  species  react  to  changing  depth  in  different  ways  and  that 
different  experimental  methods  can  affect  the  way  a fish  would  normally  react  Hence,  there  are 
uncertainties  in  the  actual  values  of  a at  depth.  Since  the  swimbladder  is  attached  to  various  other  parts  of 
the  fish,  there  is  no  reason  to  expect  a change  in  swimbladder  volume  to  cause  a uniform  change  in  its 
linear  dimensions.  Thus,  the  uncertainty  in  a at  depth  produces  a smaller  uncertainty  in  e.  Various  depth 
variations  in  a were  examined  when  comparing  the  data  to  the  model.  However,  since  the  individual 
researcher  is  in  the  best  position  to  interpret  his  own  data,  the  final  comparisons  utilized  the  depth  variations 
chosen  by  those  researchers. 

The  new  model  requires  values  of  E and  s.  Since  the  exact  value,  or  range  of  values,  of  E is  presently 
unknown  and  s is  quite  likely  under  the  control  of  the  individual  fish,  direct  comparisons  of  the  model  and 


experimental  data  are  not  possible.  However,  indirect  comparisons  are  possible.  The  first  step  in  a 
comparison  is  to  determine  the  experimental  values  of  Q,  a,  and  e.  Then  is  set  equal  to  u)b,„  and  s 
calculated  utilizing  equations  IV- 1 and  IV- 18  and  figure  1 . If  surface  tension  has  a negligible  effect  on  s 
can  be  assumed  to  be  equal  to  10*  dyne/cm,  based  on  the  data  given  in  Appendix  A.  Q,„  and  0,^  are  then 
calculated  from  equations  IV-4  and  IV-6  and  Q„  is  determined  from  the  equation; 

u " 5;:^  ^ i ^ ■ (iv-19) 

Finally,  ^ is  calculated  from  equation  IV-5.  The  results  of  this  procedure  are  given  in  Table  II,  where  L is  the 
fish  length,  D,  a,  <%„,  and  Q are  experimental  data,  and  s and  ^ are  calculated  by  equating  the  model  to  the 
data 

The  new  model  can  be  indirectly  compared  to  the  experimental  data  by  examining  the  values  of  s and  ^ 
given  in  Table  II.  The  values  shown  in  Table  II  are  within  the  limits  chosen  in  Appendix  A to  give  reasonable 
ranges  for  s and  This  is  a necessary  condition  for  the  model  to  be  valid,  but  it  is  by  no  means  a sufficient 
condition  In  order  that  the  model  can  be  used  with  some  degree  of  confidence  to  predict  resonant 
scattenng  from  swimbladdcr-bearing  fish,  the  variations  in  s and  ^ shown  in  Table  II  must  also  be  explained. 

The  five  sets  of  data  can  be  separated  into  two  groups,  based  upon  the  magnitude  of  tension  in  the 
swimbladder  wall  Swimbladder  tension  had  little  or  rxj  effect  on  the  measurements  made  by  Coate,  Batzler  and 
Pickwell.  arxl  Sundries  and  Sand,  where  s < 10^  dyne/cm.  Swimbladder  tension  appears  to  have  had  a significant 
effect  on  the  measurements  by  McCartney  and  stubbs  and  Sand  and  Hawkins,  where 
2 X 10^  dyne/cm  < s < 4 x io^  dyne/cm. 


Table  II  — Expenmenta!  Data  and  Results  of  Comparisons  to  the  New  Model 


Source 

Fish 

L 

cm 

D 

m 

a 

cm 

rad/sec 

Q 

s 

dyne/cm 

poise 

Coate  [36] 

Crappie 

20 

4.6 

1 10 

2400 

4.6 

10* 

300 

Batzler  and 

Goldfish 

7 

6 

0.51 

5650 

5,1 

6 xio* 

130 

Pickwell  [15] 

Goldfish 

6 

6 

0 44 

6600 

3.8 

5 xio* 

160 

Anchovy 

11 

6 

0.40 

8000 

4.5 

7 xio* 

130 

McCartney  and 

Coalfish 

30 

10 

1,17 

4150 

1.1 

1.7x10® 

2600 

Stubbs  [37] 

20 

1.08 

5200 

1.7 

1.9x10® 

1750 

30 

1.03 

6000 

2.5 

1.7x10® 

1200 

Pollack 

35 

30 

1.13 

4800 

1.4 

3.6x10® 

2200 

40 

1.09 

5800 

1.3 

6.1x10® 

2600 

Ling 

55 

30 

1.60 

3150 

2.5 

8 xio® 

1500 

Cod 

35 

30 

1.52 

3500 

2.0 

1.6x10® 

2000 

Sand  and 

Cod 

16 

11 

0.69 

10700 

1.0 

3.8x10® 

2500 

Hawkins  [33] 

25 

0.58 

8000 

1.3 

2.8x10® 

1000 

30 

0 55 

8800 

1.6 

2.1  xio® 

830 

35 

0.53 

9900 

2.2 

3.3x10® 

610 

40 

0,51 

10900 

2.5 

3.8x10® 

540 

45 

0.50 

12000 

3.2 

5.1  xio® 

430 

50 

0.48 

12800 

3.5 

4.7x10® 

370 

Sundries  and 

Charr 

40 

2 

0.95 

3000 

5.2 

10* 

240 

Sand  [38] 

[Averages 

6 

0.87 

4000 

4.5 

10* 

300 

for  5 

10 

0.80 

4900 

5.0 

10* 

280 

fish) 

15 

0.74 

5800 

6.1 

10* 

220 

29 


1 


aa 


Goldfish,  anchovies,  and  charr  are  physostomes.  that  is,  their  swimbladders  have  an  opening  into  their 
stomachs.  Hence,  based  on  consideration  of  equation  11-115,  it  is  to  be  expected  that  the  tension  in 
physostome  swimbladders  would  be  relatively  low.  The  data  show  that  this  is  so.  The  other  fish  are 
physoclists,  that  is,  their  swimbladders  are  completely  closed.  These  fish  may  well  obtain  some  benefit,  either 
hydrostatic  or  otherwise,  by  maintaining  a tension  in  the  swimbladder  wall  The  data  indicate  that  cod,  ling, 
pollack,  and  coalfish  maintained  relatively  high  values  of  s Physiologically,  it  should  be  possible  for  a 
physoclist  to  vary  the  tension  in  the  swimbladder  wall  from  very  taut  to  flaccid  It  appears  that  the  crappie 
swimbladder  was  in  a flaccid  condition. 

Since  physoclists  probably  can  vary  the  tension  in  their  swimbladders.  any  shock  to  their  systems,  such  as 
a rapid  increase  in  pressure,  might  cause  the  swimbladder  tension  to  change  dramatically.  This  seems  to  be 
the  situation  for  the  cod  examined  by  Sand  and  Hawkins  This  fish  had  been  allowed  to  become  adapted  to  a 
depth  of  1 1 m for  at  least  48  hours.  At  this  depth,  its  measured  resonant  frequency  was  several  times  higher 
than  that  of  a bubble  of  the  same  presumed  size  and  shape  This  difference  was  attributed  to  a high 
swimbladder  tension.  Rapidly  shifting  this  fish  to  a depth  of  25  m significantly  decreased  the  measured  resonant 
frequency  to  a value  only  about  10  percent  greater  than  that  expected  for  a bubble  of  the  same  presumed  size  and 
shape.  This  ratio  remained  essentially  constant  as  the  fish  was  shifted  in  5 m increments  to  50  m.  The 
decrease  in  as  the  fish  was  shifted  from  11  m to  25  m is  readily  explained  by  the  new  model  by  a 
sudden  decrease  in  s,  which  seems  quite  feasible  from  a physiological  standpoint 

Thus,  the  variations  in  the  values  of  s calculated  using  the  new  model  can  be  explained  by  the 
differences  between  physostomes  and  physoclists  and  by  the  ability  of  physoclists  to  vary  swimbladder 
tension.  Hence,  the  new  model  is  an  improvement  over  the  models  of  Andreeva  and  Lebedeva.  By  using  a 
variable  swimbladder  tension,  rather  than  a fixed  shear  modulus,  to  model  tissue  stiffness,  the  new  model 
can  predict  resonant  frequencies  significantly  higher  than  those  expected  for  a free  bubble,  whereas  the 
other  models  cannot 

As  a practical  matter,  a parameter  that  can  be  randomly  varied  over  several  orders  of  magnitude  is  not 
very  useful  in  a predictive  model.  However,  the  variations  in  tension  indicated  in  Table  II  were  obtained  from 
fish  which  were  subjected  to  other  than  natural  conditions,  such  as  rapid  changes  in  depth.  It  is  quite 
possible  that,  under  more  natural  conditions,  values  of  swimbladder  tension  would  be  much  more  uniform. 
Thus,  if  coukJ  very  well  be  that  further  experiments,  in  wnich  the  experimental  conditions  more  closely 
approximate  natural  conditions,  could  provide  much  information  about  swimbladder  tension.  This  would 
enhance  the  value  of  the  new  model  as  a predictor  of  resonance  in  large  physoclists 

Although  the  variations  in  s can  be  adequately  explained  by  variations  in  fish  physiology,  the  uncertain 
variation  of  a with  depth  definitely  causes  some  degree  of  variation  in  s In  order  to  calculate  s,  the  difference 
between  oib^and  oig  for  a bubble  of  the  same  assumed  size  and  shape  is  required  Thus,  the  accuracy  of  the 
estimate  of  a greatly  affects  the  accuracy  of  the  calculated  value  of  s 

Sundnes  and  Sand  measured  the  resonant  frequencies  of  charr  in  order  to  determine  their  swimbladder 
volumes.  In  the  case  of  these  physostomes,  this  was  a valid  procedure  However,  in  the  case  of  large 
physoclists,  where  swimbladder  tension  affects  a>n,„,  this  procedure  would  be  invalid  Hence,  experiments 
which  might  be  designated  to  acoustically  examine  swimbladder  tension  require  an  accurate  and 
independent  measurement  of  swimbladder  volume  at  depth 

The  values  of  ^ given  in  Table  II  are  also  grouped  on  the  basis  of  the  magnitude  of  swimbladder  tension. 
For  the  cases  where 

s < 10*  dyne/cm, 

the  values  of  ^ range  from  130  to  300  poise,  with  the  values  being  quite  consistent  for  a particular  set  of  data 
For  the  cases  where 

s > 2 X 10*  dyne/cm, 

the  values  of  E,  range  from  370  to  2,600  poise,  with  wide  variations  within  a particular  set  of  data  Based  on  the 
consistent  results  obtained,  and  the  reasonableness  of  their  values,  it  appears  that  the  new  model  can  be 
used  to  predict  the  damping  for  the  first  group.  However,  the  high  values  and  variability  of  ^ in  the  second 
group  must  be  examined  more  closely  before  any  conclusion  can  be  reached 
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There  are  essentially  two  sources  of  the  variations  in  the  calculated  values  of  ^ between  the  two  groups 
and  within  the  group  where  s > 2 x io^  dyne/cm.  One  source  is  the  experimental  data  and  the  other  is  the 
model.  The  most  likely,  and  probably  greatest,  source  of  error  in  the  experimental  data  is  the  uncertainty  in 
swimbladder  volume.  Since  the  calculated  values  of  ^ are  proportional  to  a^,  an  error  in  a can  cause  a large 
error  in  However,  as  mentioned  earlier,  the  experimental  data  were  examined  assuming  various 
swimbladder  volume-depth  relationships.  None  of  these  relationships  produced  consistent  results  or  values 
as  low  as  those  obtained  for  the  first  group,  although  in  some  cases  the  variability  was  reduced.  For  example, 
if  a constant  a=0.69  cm  is  assumed  for  the  cod  examined  by  Sand  and  Hawkins,  ^ varies  from  1 ,500  to  900 
poise  as  the  depth  increases  from  25  to  50  m.  This  variation  is  smaller  than  that  for  the  swimbladder-depth 
relationship  which  was  utilized,  but  the  values  are  still  much  higher  than  those  of  the  first  group  and  are  by  no 
means  consistent.  It  is  possible,  though,  that  the  cod-like  fish,  which  have  a well-developed  muscle  system, 
have  an  intrinsically  higher  tissue  viscosity  than  fish  such  as  goldfish  or  anchovies. 

At  the  beginning  of  this  chapter,  it  was  assumed  that 
^<2x10^  poise 

and  the  equations  were  simplified  accordingly.  It  might  be  suspected  that  this  could  be  a source  of  error  in  the 
cases  where  ^ > 2 x W poise.  However,  the  error  caused  by  simplifying  equations  III-26  and  III-34  to 
equations  IV- 1 and  IV-2  is  less  than  ten  percent  for  values  as  high  as  ^ = 6 x ICP  poise.  Thus,  the  simplified 
equations  are  not  a significant  source  of  error. 

The  other  primary  source  of  error  is  the  model  itself.  Since  it  does  not  seem  to  be  possible  to  attribute  all 
the  variations  in  ^ to  the  errors  in  the  data,  it  appears  that  the  remaining  variations  are  due  to  inadequacies 
in  the  model. 

An  examination  of  all  the  data  in  Table  II  indicates  that  ^ generally  increases  with  s.  There  is  no  indication 
that  this  is  so  for  the  group  of  data  for  which  s < 1 0^dyne/cm.  However,  there  is  definite  correlation  between  s 
and  ^ when  the  two  groups  of  data  are  compared.  Also,  despite  the  wide  variations  of  ^ within  the  group  of 
data  for  which  s > 2 x tCF  dyne/cm,  there  is  a rough  indication  that  E,  increases  with  s for  this  group  of  data. 
Thus,  it  appears  that  increasing  tension  in  the  swimbladder  wall  has  little  or  no  effect  on  ^ for  s < 1 0^  dyne/ cm, 
but  that  it  does  cause  E,  to  increase  tor  s > 2 x 10^  dyne/cm.  The  apparent  increase  in  viscosity  with  increasing 
tension  in  the  swimbladder  wall  for  s > 2 x 1 0^dyne/cm  would  not  be  expected  to  occur  in  a Newtonian 
fluid.  This  indicates  that  modelling  a fish  as  a Newtonian  fluid  may  not  be  appropriate  for  s > 2 x 10'' 
dyne/cm  and  that  some  type  of  non-Newtonian  model,  such  as  a dilatant  or  viscoelastic  fluid,  may  be  more 
appropriate.  However,  the  present  model  is  appropriate  tor  s <10^  dyne/cm. 
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CHAPTER  V 


CONCLUSIONS 


The  new  model  of  a swimbladder-bearing  fish  as  a viscous  shell  enclosing  an  air  cavity  with  surface 
tension  at  the  inner  interface  was  developed  because  older  models,  of  which  Andreeva's  model  [26]  is  most 
widely  employed,  do  not  accurately  predict  the  results  obtained  in  experiments  on  the  characteristics  of 
swimbladder  resonance.  The  principal  difference  between  the  older  models  and  the  data  is  that  the  values 
of  Q predicted  by  the  models  are  always  higher  than  those  obtained  experimentally.  A second  difference  is 
that  the  older  models  can  not  account  for  the  high  resonant  frequencies  obtained  for  some  large 
physoclists.  The  new  model  sought  to  correct  these  differences  by  explicitly  including  the  viscosity  of  fish 
tissue  and  by  including  a tension  in  the  swimbladder  wall. 


For  convenience,  the  equations  developed  for  the  new  model  will  be  restated  here: 
a,o^a=  = ^ ^ (3Y.  - 1)  . 
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Comparable  equations  for  an  air  bubble  in  water  are  available  at  only  co  = cOq,  where  H = Q.  The  equations 
for  o and  Q obtained  for  the  new  model  are  essentially  the  same  as  those  for  an  air  bubble  in  water,  the  only 
significant  difference  being  in  the  value  of  the  viscosities  of  fish  flesh  and  water.  However,  the  equations  for 
CJti  for  the  new  model  and  an  air  bubble  in  water  for  which  viscosity  and  surface  tension  are  included  are  not 
the  same.  If  the  viscosity  of  fish  flesh  were  used  in  the  equation  obtained  for  an  air  bubble,  ojo  could  be  zero 
or  imaginary.  However,  for  the  new  model,  the  effect  of  viscosity  on  cjo  was  shown  to  be  small  enough  to  be 
neglected. 

Equations  V-1  and  V-2  are  valid  for  an  upper  limit  for  ^ = 6 x ICF  poise.  If  ^ > 6 x 10^  poise,  then  the 
outer  shell  radius,  b,  appears  in  the  o-equation.  This  implies  that  this  is  a boundary  layer-type  problem. 
Thus,  if  the  fish  flesh-water  interface  is  outside  the  boundary  layer,  the  magnitude  of  b is  immaterial. 

All  of  the  models  are  spherical  in  nature.  However,  the  swimbladders  of  some  fish  are  sufficiently 
elongated  that  their  shapes  can  have  a significant  effect  on  cOq.  Thus,  equation  V-1  should  be  modified  to 
include  this  effect.  Hence: 
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Although  Z,  was  determined  for  a bubble  in  water,  its  use  here  is  quite  reasonable,  especially  when  the 
similarity  between  the  equations  for  a bubble  and  the  new  model  are  considered. 

Several  other  conclusions  can  be  reached  for  the  new  model.  One  is  that,  for  the  ratios  of  outside  to 
inside  shell  diameters  considered,  the  actual  value  of  the  outside  diameter  has  no  effect  on  the  results. 
Another  is  that  thermal  losses  are  not  very  significant  and,  in  most  cases,  can  be  neglected.  A third  is  that, 
for  low  values  of  H,  ((Oq/wH)^  can  be  comparable  to  [(cDo^/o)^)  - If  at  off-resonant  frequencies.  Hence, 
considering  the  experimental  values  of  Q obtained,  it  is  apparent  that  H,  rather  than  Q,  should  be  used  to 
calculate  o. 

A comparison  of  the  new  model  with  available  experimental  data  indicates  that  the  new  model 
constitutes  a definite  improvement  over  previous  models.  The  new  model  can  predict  low  values  of  Q and 
elevated  values  of  tOo  which  the  previous  models  could  not.  In  addition,  the  new  model  can  be  used  to 
obtain  the  magnitude  of  damping  at  any  frequency,  whereas  many  previous  models  only  produced  the 
value  at  resonance. 
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The  comparison  of  the  new  model  with  experimental  data  indicates  that  the  model  is  most  accurate  tor 
fish  in  which  s < 1 0’ dyne/cm.  Table  II  indicates  that  this  includes  physostomes  up  to  at  least  40  cm  in 
length.  In  addition,  Appendix  A indicates  that  small  physoclists.  10  cm  or  less  in  length,  are  also  included 
in  this  group.  Most  of  the  fish  found  in  deep  scattering  layers  are  smaller  than  1 0 cm  and  none  have  the  well 
developed  musculature  of  fish  such  as  cod.  Thus,  it  appears  that  the  new  model  will  be  of  considerable 
value  in  studies  of  volume  reverberation.  Table  II  indicates  that  a value  of  ^ of  approximately  200  poise  is 
appropriate  for  these  studies 

The  variations  in  s and  E,  reguired  to  match  the  values  calculated  using  the  new  model  to  the 
experimental  data  for  large  physoclists  indicate  that  the  new  model  is  not  completely  adequate  for  these 
fish.  Even  so,  the  new  model  is  still  better  than  previous  models  and,  as  such,  is  of  some  value  in  studies  of 
resonance  of  targe  physoclists.  Table  II  indicates  that  values  of  s = 1 0’  dyne/cm  and  1 0’  poise  5^22x10’ 
poise  should  give  reasonable  results  for  these  studies. 

The  new  model  introduces  two  new  parameters,  E,  and  s.  Few  measurements  of  these  parameters  exist 
and  those  that  do  vary  widely.  The  accuracy  of  the  new  model  could  be  determined  with  much  more 
confidence  if  better  information  on  ^ and  s were  available.  In  addition  the  variations  of  swimbladder  volume 
with  depth  for  different  fish  species  are  required  for  all  models.  Hence,  further  experiments  to  determine  s 
and  swimbladder  size  and  shape  versus  depth  are  recommended.  This  recommendation  is  made  with 
the  realization  that  the  complexities  involved  in  these  experiments  will  be  significant. 
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APPENDIX  A 

PHYSICAL  PROPERTIES  OF  FISH 


Several  of  the  physical  properties  of  fish  flesh  which  must  be  known  for  this  study  have  not  been 
accurately  measured.  In  this  appendix,  the  available  data  will  be  discussed  and,  where  necessary, 
appropriate  approximations  determined.  Properties  which  must  be  specified  include  compressional  wave 
velocity  (or  sound  velocity),  c,;  density,  po,;  specific  heat  at  constant  pressure,  Cp,;  ratio  of  specific  heats,  Vii 
thermal  conductivity,  shear  viscosity,  ri^;  and  bulk  viscosity,  r|  . In  addition,  the  ratio  of  swimbladder 
volume  to  total  fish  volume  and  the  surface  tension,  s,  at  the  air-fish  flesh  interface  must  be  determined. 

Experimental  evidence  indicates  that  swimbladder  volumes  of  small  mid-water  fishes  range  from  about 
0.5  to  5 percent  of  the  total  fish  volume  [32].  For  larger,  near  surface  marine  fishes,  swimbladder  volumes 
are  about  4 to  5 percent  of  the  total  fish  volume  [At].  For  small  fish,  the  ratio  of  outer  radius,  b,  to  inner 
radius,  a,  of  the  fish  flesh  shell  is  chosen  to  be  2.5  Sb/a  S6,  which  corresponds  to  swimbladder  percentage 
volumes  of  about  6 to  0,5  percent.  For  large  fish,  b/a  is  chosen  to  be  2.5  <b/a  <3.2,  which  corresponds  to 
swimbladder  percentage  volumes  of  about  6 to  3 percent.  The  fish  size  range  of  interest  is  about  1 cm  to  1 m, 
which  roughly  corresponds  to  0. 1 cm  < a < 5 cm  [32.  A2,  A3]. 

The  acoustic  properties  of  fish  flesh  have  been  measured  by  several  researchers  [31,  A1.  A4-A61. 
Experimental  values  of  density  range  from  about  1 .02  to  1 .09  gm/cm^,  with  an  average  of  about  1 .05 
gm/cmT  Experimental  sound  velocities  range  from  about  1.50  x lO^  to  1.60  x 10^  cm/sec,  with  an  average 
of  1.55  X 10®  cm/sec.  The  average  values  are  used  in  this  report. 

Measurements  on  the  thermal  properties  of  fish  flesh  are  not  available.  However,  data  on  the  thermal 
conductivity  and  specific  heat  at  constant  pressure  of  human  and  dog  tissue  do  exist  [A-7].  These  data 
provide  sufficiently  close  approximations  to  the  required  values.  Thus,  for  fish  tissue 


c 


PI 


0.89 


cal 

gm°C 


and 


cal 


K,  = 1.32  X 10  - „„ 

’ cm  sec  C 

These  values  are  quite  close  to  those  for  sea  water  of  35  parts  per  thousand  salinity,  which  are  [45] 


cal 


and 


Cp.  - 0.93 

= 1.34  X 10  3 


cal 


Also  for  sea  water  [45], 


cm  sec  °C 


Y.  » 1.01  . 

Thus,  it  will  be  assumed  that,  since  other  thermal  properties  of  flesh  and  sea  water  are  very  similar, 


Y,  « 1.01  . 


Both  the  shear  viscosity,  ris„  and  the  bulk  viscosity,  ri^,,  are  required  for  the  present  model.  However,  for 
convenience,  a viscosity  parameter,  L will  be  defined  as 

^ = 4/3  ri3, -I- n,, . (A-1) 

The  ratio,  Pf/ris,  for  animal  tissue  is  similar  to  that  for  water,  which  is  approximately  3 [A8.  A9)  Thus, 

^ 2:4.3r|„  s?  1.4  ri^,  . ^^.2) 

Only  one  set  of  data  on  the  viscosity  of  animal  tissue  is  available.  However,  other  data  exist  from  which 
tissue  viscosity  can  be  determined  indirectly.  These  data  are  measurements  of  absorption,  complex  shear 
modulus,  and  cell  viscosity  of  animal  tissue.  In  all  cases  it  will  be  assumed  that  the  viscosity  of  all  animal 
flesh  is  approximately  equal. 

The  direct  measurements  of  viscosity  v'pre  performed  on  mammalian  tissue  [A  10]  h,,  was  determined 
by  four  different  methods  The  results  ranged  from  100  to  420  poise,  with  an  average  value  of  175  poise 
Thus  the  values  for  ^ range  from  430  to  1800  poise,  with  an  average  value  of  760  poise. 


Viscosity  can  be  determined  from  measurements  of  absorption  of  sound  by  utilizing  the  equation  [All]; 

F = (A-3) 

^ 0)2 

where  a is  the  absorption  and  o)  the  circular  frequency  Many  measurements  of  absorption  in  tissue  have 
been  made  at  high  frequencies.  However,  tissue  exhibits  relaxation  phenomena  at  high  frequencies  so  that 
viscosities  at  frequencies  above  relaxation  cannot  be  directly  related  to  those  below  and  can  be  used  only 
as  lower  limits.  It  has  been  found  that  for  muscle,  q,  relaxes  near  400  kHz  and  r|g  relaxes  near  40  kHz 
(A  12].  Only  one  set  of  absorption  measurements  has  been  conducted  below  the  Megahertz  range.  These 
measurements  were  made  at  300  to  350  kHz  [A  13).  Utilizing  equation  A-3,  ^ was  calculated  from  these  data 
to  range  from  30  to  220  poise,  with  an  average  value  of  130  poise.  However,  ri^  has  already  relaxed  at  the 
measurement  frequencies,  so  that  its  contribution  to  the  absorption  is  unknown  Thus,  equation  A-3 
provides  a lower  limit  to  An  upper  limit  can  be  determined  for  E,  if  it  is  assumed  that  the  absorption  is  due 
solely  to  ris.  Then,  at  lower  frequencies,  ^ would  range  from  1 30  to  950  poise,  with  an  average  value  of  580 
poise.  Thus,  from  the  absorption  measurements: 

30  poise  < ^ < 950  poise. 

Viscosity  can  also  be  estimated  from  measurements  of  complex  shear  modulus  by  utilizing  the  equation 
[A9]: 

Is  =-^  . (A-4) 

where  p,  is  the  imaginary  part  of  the  complex  shear  modulus.  One  set  of  measurements  of  the  complex 
shear  modulus  of  fish  tissue  has  been  made  from  about  2 to  14  kHz  [27].  Utilizing  equations  A-4  and  A-1 . ^ 
was  calculated  from  these  data  to  range  from  about  4 to  90  poise,  with  an  average  value  of  28  poise. 

The  viscosity  of  animal  tissue  can  also  be  estimated  from  the  viscosity  of  animal  cells  if  it  is  assumed 
that  the  tissue  viscosity  is  equivalent  to  the  cell  viscosity.  Separate  measurements  have  been  made  on  the 
viscosities  of  both  cell  protoplasm  and  membrane.  Thus,  to  estimate  the  viscosity  of  the  complete  cell,  a 
geometric  average  of  the  membrane  and  protoplasmic  viscosities  is  calculated  based  on  the  proportional 
thicknesses  of  membrane  and  protoplasm.  The  equation  used  to  calculate  cell  viscosity  is 

F = ra(b  — a)  ^ c 1 ’ 

^ [ b2  J ’ (A-5) 

where  E^,  ^ and  ^ are  the  viscosity  parameters  of  the  cell,  protoplasm  and  membrane,  respectively,  and  a 
and  b are  the  inner  and  outer  membrane  radii.  Cell  radii  range  from  2 x 10  ■*  to  15  x 10  “ cm  and  cell 
membranes  are  75  x 10  ® to  10  ® cm  thick  [A  14],  Thus 

5 X 10  ^ <-5^  <5  X 10  3, 
b 

and  a/b~  1 . Measurements  have  been  made  on  ris  of  protoplasm  and  r^p  of  membranes.  For  protoplasm,  ri^ 
was  found  to  range  from  4 x 10  2 to  3 x 10  ' poise  [A15].  For  membranes,  was  found  to  range  from  2.7 
X IOM0  2.7  X 10®  poise  [A16].  Therefore,  utilizing  equations  A-2  and  A-5, 

60  poise  < - 1.600  poise. 

This  estimate  is  probably  subject  to  the  greatest  error  of  the  three  indirect  estimates  of  viscosity  due  to  all 
the  assumptions  required. 

Summarizing  the  ranges  of  E determined  by  the  various  methods  in  their  probable  order  of  accuracy; 
direct  measurement, 

430  poise  5 ^ S 1,800  poise; 
absorption, 

30  poise  < ^ S 950  poise; 
complex  shear  modulus, 

4 poise  < ^ < 90  poise; 
cell  viscosity, 

60  poise  S ^ < 1,600  poise. 


The  total  range  of  data  is 

4 poise  < ^ < 1 800  poise 
Thus,  for  this  study,  limits  of 

1 poise  < 5 < 10“  poise 
will  be  used,  with  a more  likely  range  being 
50  poise  <^<2x10^  poise. 

The  swimbladder  wall  of  a fish  is  a membrane  which  is  capable  of  supporting  tension.  In  the  present 
model,  the  swimbladder  wall  has  zero  thickness,  so  that  the  tension  in  the  membrane  is  effectively  a surface 
tension.  Measurements  of  the  internal  swimbladder  pressure  can  be  used  to  calculate  surface  tension 
since 

AP  = (2s/a),  (A-6) 

where  AP  is  the  difference  between  the  internal  swimbladder  pressure  and  the  ambient  pressure  [50]. 
Several  researchers  have  measured  internal  swimbladder  pressures,  but  since  it  is  probable  that  a fish  can 
control  the  tension  in  the  swimbladder  wall,  only  measurements  on  live,  unanesthetized,  fish  are 
considered.  Excess  internal  pressures  from  2 x 10-'  to  6 x io^  dynes/cm^  have  been  measured  in  fishes 
which  were  2 to  20  cm  long  [Al?].  These  results  correspond  to  surface  tensions  of  6 x lO'*  to  7 x 10“ 
dyne/cm  for  fish  with  swimbladder  radii  of  about  0.1  to  1.0  cm.  The  majority  of  these  fish  were  less  than  10 
cm  long  with  swimbladder  radii  less  than  0.5  cm. 

Surface  tension  is  probably  a function  of  fish  size,  so  that  measurements  on  larger  fish  are  needed  for 
the  present  study.  No  measurements  of  excess  pressure  of  larger,  unanesthetized  fish  are  available. 
However,  another  means  can  be  used  to  estimate  the  upper  limits  of  surface  tension  in  larger  fishes.  As 
discussed  in  the  text.  Sand  and  Hawkins  have  attributed  high  experimental  resonant  frequencies  to 
swimbladder  tension  [33].  If  this  is  true,  then  an  upper  limit  of  surface  tension  can  be  calculated  by 
assuming  that  the  fish  is  a free  bubble  with  surface  tension  and  utilizing  equation  1-12,  neglecting  viscosity, 
and  equation  1-25  to  account  for  spheriodal  swimbladder  shapes.  Although  this  method  does  not 
necessarily  give  accurate  estimates  of  surface  tension,  it  does  produce  values  which  can  be  used  as  upper 
limits.  These  limits  are  useful  because  it  is  the  possible  range  of  surface  tension  that  is  required.  Surface 
tensions  calculated  from  resonance  measurements  range  from  about  ICT  to  10®  dyne/cm  for  swimbladder 
radii  from  1 to  2.5  cm  [33,37]. 

[ i The  surface  tension  of  an  air  bubble  in  water  is  74  dyne/cm  [50].  Hence,  the  range  of  surface  tension  for 

f ‘ small  swimbladders  (aw  0.1  cm;  is  chosen  to  be 

|.  10^  dyne/cm  5 s S 10®  dyne/cm. 

r 

I For  larger  swimbladders  (a si  5 cm)  the  range  is  chosen  to  be 

I 10^  dyne/cm  - s ^ 10®  dyne/cm. 

I ' 
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APPENDIX  B 

SIMPLIFIED  EXPRESSION  FOR  (Bw/A) 


Simplification  of  (B„/A)  as  given  in  equation  111-16  is  based  on  the  ranges  of  physical  properties  and  the 
limits  of  the  variables  given  in  Chapter  11. 

The  first  step  in  the  simplification  of  equation  111-16  is  the  simplification  of  the  spherical  Bessel. 
Neumann,  and  Hankel  functions  in  the  S,|’s.  In  the  ki^a  terms: 


(B-1) 

and 

. ,,,,  cos(k,,a)  sin(k,„a) 

(B-2) 

Now,  if 

sin(k,,a)  = sin[(  1 -i-i)  xa]  , 

(B-3) 

then 

sin(k,,a)  = sin(xa)cosh(xa)  - i sinh(xa)cos(xa) 

(B-4) 

Similarly, 

cos(k,3a)  = cos(xa)cosh(xa)  - i sin(xa)sinh(xa)  . 

(B-5) 

An  examination  of  the  parameters  involved  shows  that  xa  > 1 0,  so  that. 

cosh(xa)  sinh(xa)  ;;; 

(B-6) 

and 

sin(k,3a)  s:  i cos(k,,a)  . 

(B-7) 

Thus, 

Wk,.a)  - 

(B-8) 

and 

la” 

jo'(k,3a)  = -jo(k,aa)  [i  + ^ ] . 

(B-9) 

One  of  the  assumptions  in  the  model  is  that  the  shell  is  small  compared  to  the  wavelength  of  the  incident 
compressional  wave.  This  means  that  kj^b,  k^jb,  and  kj^a  are  small.  The  definition  of  "small " will  now  be 
determined  by  examining  the  expansions  for  spherical  functions  of  small  argument: 


jo(z)  = 

(’- 

). 

(B-10) 

jo'(z)  = 

z 

3 

('-4--  )■ 

(B-11) 

jo"(z)  = 

-4' 

(B-1 2) 

no(z)  = 

1 1 

(B-1 3) 

-t' 

[1  --y  + . . . . 

no'(z)  = 

1 1 

4 Z2  \ 

(B-1 4) 

Z2  i 

1 +“2  - • • ■ • 

no"(z)  = 

2 

Z3 

1 

-1- 

(B-1 5) 
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- X 


A 


• « 
-A 


ho(z)  = I 

Z2 

6 

(6-16) 

ho'(z)  = 

z 

T 

0 

Z2 

-To  ^ 

)■ 

(B-17) 

ho''(z)  = 

1 

3 

(' 

10 

■ ) 

(B-18) 

In  order  to  limit  errors  in  these  functions  to  under  ten  percent,  it  must  be  assumed  that 


< 2x10'  , (B-19) 

k2,2b^<  2 X 10  ' (B-20) 

and 

kj^^a^  < 6 X 10  ’ . (B-21) 

This  implies  that  in  fish  flesh  and  water 

cob  < 7 X 10''  cm/sec  (B-22) 

and  in  air 

coa  < 2.4  X 10''  cm/sec  . (B-23) 


Then  with  these  limitations  and  the  acceptance  of  at  most  a ten  percent  error,  equations  B- 10  through  B- 18 
can  be  written  as: 


io(z)  = 1 , 

{B-24) 

)o'(Z)  = 

(B-25) 

jo"(z)  = -y  . 

(B-26) 

no(z)  = - — , 

(B-27) 

no'{z)=^  , 

(B-28) 

no"(z)  = -y  . 

(B-29) 

ho(z)  = 1 -y  . 

(B-30) 

ho'(z)  = + ^ . 

{B-31) 

h„"(z)  = -l-|  . (B-32) 

It  is  necessary  to  keep  both  the  largest  real  and  largest  imaginary  term  in  the  spherical  Hankel  function 
expansions,  regardless  of  their  comparative  magnitudes. 

Equations  B-24  through  B-32  will  now  be  utilized  to  simplify  equations  III-2  through  III-15.  In  addition, 
equations  11-59, 11-60. 11-68,  and  11-69,  will  be  substituted  for  k,,.  kj,,  k^,,  and  kj,,.  These  substitutions  yield, 
after  neglecting  any  term  which  causes  an  error  of  less  than  ten  percent,  the  following  expressions  for  the 
S.;s: 

£,  w2b  ic* 

” PowWb2 

(B-33) 
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where 


fa)^b 

3po»c,2 


(B-49) 


pCb)^-"] 


(B-50) 


The  remainder  of  the  simplification  process  for  equation  111-16  will  be  conducted  in  several  steps.  In  each 
step,  the  order  of  magnitude  of  the  various  terms  will  be  compared  and  terms  which  cause  an  error  of  less  than 
ten  percent  will  be  neglected.  This  procedure  will  be  done  in  such  a way  that  no  term  is  neglected  which  could 
become  important  due  to  a later  subtraction.  It  should  be  mentioned  that  most  of  the  terms  which  are 
neglected  result  in  errors  of  much  less  than  ten  percent  and  in  those  cases  where  the  errors  approach  ten 
percent,  that  error  is  approached  only  at  the  limiting  values  of  the  terms  involved.  In  order  that  the 
simplification  process  can  be  checked  by  the  interested  reader,  the  results  of  each  step  will  be  given,  rather 
than  just  the  final  result. 

The  next  step  is  to  simplify  the  terms  in  the  parentheses  in  equations  111-17  and  1(1-18.  The  results  are 


S22S33  823832  - ( ) ( 


6Po,^c,aJ 


(B-51) 


S23S42  S22843  — 


(B-52) 


8,3832  8,2833  - ( 3p 2 ) (1  ir)-  ' (2po,^S=  ) (^  -fe- ) 


(B-53) 


8,2843  8,3842  - (^p^)  + ( 

■ (p;;i3br)  ^(epc^lb^ ) ) 


(B-54) 


The  next  apparent  step  would  be  to  simplify  the  expressions  for  U and  W.  This,  however,  can  lead  to  errors 
when  the  subtractions  indicated  for  Y and  Z are  made.  Hence,  the  next  step  is  to  simplify  Y and  (Y-Z).  Before 
writing  the  simplified  expressions  for  Y and  (Y-Z),  several  intermediate  steps  will  be  given  in  order  to  assist  the 
interested  reader  in  checking  the  final  expressions. 
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Substitution  of  equations  B-43,  B-46,  B-51,  and  B-52  into  III-17  yields: 

/ c,  \ . 2s  \ r 2sc,(y.-1)  1 / (DK. 


\ Po,wa^)  \ 3po,c/a  T I PcPo.c.^oi^a*  J ( 2Po.Cp,  ) 

* (p^T-a)  l(^)(’-i)‘(^)(’-|^) 


- ’•  (^)"(’  -i)  - [W]  (^)“(' -w) 


rii^^(Ya-i)i 

( WK.  \Y 

l3Po, 20,232  J 

V2Po.Cpa  / V 

' ( 6po,2c,a2  ) C 3po.c.2a)  ( ^ b3  ) 

j r„  ^s(Ya-l)  -I  / ^3 

l3po,C32po,2c,(Da<  J V 2po,Cp3/  V b3  ^ 

j r-2!<YiZll!l£_l  n + /_2!_W 1 

IPo.Ca^PofCp.a^J  L M2po.u>K3^  J 

r2^(Ya-1)l  / WKa  \'^  f , 33  \ 

I Po,a)a3  J V2P03CP3/  V b3  / 

r11^^(Ya-1)1/  COK3  X''^  / a3  V I 

I 3Po,3C,3a3  J UpOaCpa  / I b3  j ) ' 


(B-55) 


Substitution  of  equations  B-43,  B-46,  B-53,  and  B-54  into  111-18  yields: 

W = (.  \ h -ii  \ ( - /'  1 - 2s  X r 2s(Ya-1)  I ( c»K,  x'^ 

V3po,3C,a2  ) V V 3Po3C.3a  ^ I PoaCa^COa^  J ^ ^PoaCp,  ) 


r 3a)^^K.(Ya-1)  1 r /2s  \ / Cp,  \ w 1 
I 2Po.C.2Po,C,2Cp3  J I Vwa^/  V 2Po.C»)Ka/  J I 

\ PoaCa'j  I V dpo,c,  ) ^ \ 3Po,3C,3a3  j V 3Po,b3  / 

r a)b(Y.-1)i/  UK.  X** 


3po,c,a  JV2Po.Cp, ) 
rU^"b(Ya-1)l  / MK,  ^ _ r C,(Ya-1)]  / “Ka  Y" 

I Po,"c,3a3  J bpo.Cp.j  I Po,wb^'  J v 2po.Cp. ) 


Po,^c,3a3 
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[ 6po,2c,b2  J I 2po,Cp.  j V as  ^ f 

K3p^)(’-&)I(4S-)(’-^) 


3^s(Y3-1)  1 / COK, 

Po.CaSpo,c,2a2  J \ 2po,Cpp 


'^•-’'(p5^)  [’■  (^)(555^)”1  1 

. / 0)  V j / (i)^a  V / 4bs  \ ra)b(Y.  - 1 ) -1  / (dk^ 

'(PoaCaV  1 V18po,^c,bs  ; i as  j I 3Po,c,a  Ji2po,Cp,  ; 

, r ^(Ya-1)  1 / X*"  /.a  , 4-  r W^^b(Ya-1)l  / WKp 

I 6Po,2C,b2  J UpppCp,  / as  ; I Po,SC,Sa3  J l2PoaCpa/ 


6Po,2c,b2  J V2P03CP 


_ r Ct(Ya-1)'j/  \ 

[ Po,U)b2  JUPoaCpa/  ) 


(B-56) 


Substitution  of  equations  B-55  and  B-56  into  B-49  yields: 


= -(^Ss;-)  (’  -^a)  ^ (^)  ( ' -f) 


M bs^  l9po,c,j  V 3po,bS;/ 


/PoaCa^Wbx  / 

2s  V 

V3po,sc,as  ) V 

3PoaCa"a; 

11^2(DSb 

27Po,Ca,2Po,2C,aS 


)(’-^)-(3ggiT)^[^]( 


Cl(Ya-1)  1 / 


2c,sb(Ya— 1)  1 / C»Ka  X r a)Sc,b(Y. - 1 ) 1 / OSK,  __a_x 

3Po„Ca,spo,a*  J V 2po,Cp,  ) [ 3po^c*2a  J \ 2poaCp,  j ^ b j 


2sb(Ya-1)l  / t«»<a  \ ^ 

3po,sc,a'  JUpoaCp,  / I bs 


_-2L\_  rm^b(Ya-1)i/  Ci)K3  \ 
^ V bs  ^ I 3po,c,a  JV2po,Cp,/ 


(»3^K,b(Y.-1) 

'l8po,c«spo,2c,Cp,a2 


I [-(^)(T55srn  (”*^) 


L 9po«c«2po,"c,a3  J V2po,Cp,  /V  b3  / 

r2a>^c,b(Y,-1)1  / \ ^ 

I 3po«c«2po,a3  J V 2Po,Cp3  /V  b3  / 

_ ra>3^K.b(Y,-1)i  r ^ /^x  / c^x'^  i _^X 

I 2po,3c,3Cp3a2  J I Vu)a2/  V2po,cuK3  ^ JV  b3  ; 


Po,3c,3a3 

PoaCa^W^^b 


1/  a)Ka  X 

Va 

I / a)K3  X 

I I 

fl1  — ^ 

1 V2PoaCpa  ) 

1 ^ 

6po,3c,b2  J 

I V 2PoaCpa  ) 

1 33  ; 

,/  Po^Ca^CD^^b  X / 2s X / ^^X 

\ 18po„c*2po,2c,a2  / V 3poaC,3a  /V  b3  / 

■ (w) 


_ i^_P2aCafW^\  /'l  —lL—\  h 

V 2po,3c,3a2  ) \ 3poaCa3a  ) \ / 


■ C,(Ya-1)l, 

( i 

r a)2b(Ya-1)l 

( WK3  X 

. Po,b^  J 

^ 2PoaCpa  / 

L 3po,c,a  J 

V2p0aCpa/ 

ra)3c,b(Ya-l)]  / V ^ 1 - — ^ 

1 3po,c«2a  J V 2po,Cp3  ) \ b ) 

[-  a)3C,Kab(Ya-'l)  1 r / 2s  X / Cp3  1 

l3po„c„2Cp3a3po,  J L Vwa3  M2PoaW'<a  ^ J 


] ['-(^)  (^)’I(’-f) 


ra)3Kab(Ya-i: 
I 3po,2c,Cp3a2 


r a)^sb(Ya-l)]  ( MKa  X . ifL') 

l9Po*c«2po,3c,a‘'J  V 2poaCp3  / V b3  / 

pa)^c,b(Ya-1)]  / tDKa  x'^  / ^ _ a3  X 
I 3po.c„3po,a3  J V 2PoaCp,  ; \ b3  ; 

r lm)3^2b(Ya-1)  1 / (iJKa  \'^ 

L 9Po»c«2po,"c.a3  J V 2poaCp3  /V  b3  ^ 


r a)^sb(Ya- 

1)1  f WKa  '1 

1 Po,3c,3a3 

J V2po,Cp3> 

IV  b' ; 

r a)^(Ya-1) 

1 ( V 

r(o2^^b(Y,-1)  1 

L epc^Qb-” 

J V2p0aCp.  f 

1"  a3  j ' 

1 Po, 30,3  as  .1 

1 ' 2P03CP3  / 

(B-57) 
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Substitution  of  equations  B-55  and  6-56  into  B-50  yields: 

PoaCa^^C. 


4a3 


/ 2^c^c,  V /.  £1\  , /_w|c^a_\  / 4b3  \ 

V 3‘*JPo„Po,azb2  ) V b3  / Vl8po,2c,b3/  1.  as  / 


( 


Po.CagCO^C^ 


, ■ V : . ) (i 


2s 


HI  + r^C,(Ya-1)]/  t4K, 

bs  / I 3Po,c,a  JV2po3CpJ 


WPoibs 
CwC,K3(Ya-1) 


L>"KOwMot^Pa 

■WCwKa(Ya-l) 


11/  tOKa  \'3 

r C»C,(Ya-1)l/ 

OJKa 

J'  2Po3Cp3  / 

1 wpo^abs  Jv 

2PoaCpa  / V 

-1)  1 Fl  + 

F 2s  w Cp3 

a^bs  J L ' 

\ was  / V2p03(DK3 

) J 

] ['  * (l^)  (t^Y  ](’  -^) 


as 


.3po,2c,Cp3as 


r 


^C«S(Ya-1) 


[ 


L3w2Po„Po,2C,a‘'b2  J V 2PoaCpa 

2^C«C,(Ya-1)  1 / CDKa 


(DK,  \ , 4as  X 

11 


bs 


[ 


“i^PowPofasbs 
11^"C„(Ya-1)l/  COK, 


](^r  (’--&■) 


3po«Po,=c,asbs  JV2poaCp3 

r^C^S(Ya— 1)1/  t4K3  / as  \ 

I Po,3c,3a^  JV2P03CP3  ^ V bs  / 

r^Ca,(Ya-1)1  ( <*>><3  4bS\ 

L6po,sc,bs  JV2P03CP3/  V ' as; 

L Po,3c,3a3  J UpoaCpa  ; V w^PowPo.a^bs  ) \ 3Po.c.sa  j 

'(■^)  (’ -4-) ('-#-) 

_ + i \ + i ( Ui!£= ) - i 

V9po,c,/  V3po,b3  ; V9powPo('c,asb3  / V bs  / V3po,3c,3as/ 


+ I 
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f i 


rc.c,(Y.-i)i 

( T - i 

r(Dc.(Y.-l)] 

1 wPo,b^  -1 

V 2po,Cp,  / 

3po,c,a  J 

V2Po,Cp.  ) 

r 2c»c,s(Ya-1)l 

( V i 

r 2c»S(Ya-1)] 

1/  UJK,  >1 

i-  (i)3po,Po«a‘'b3  J 

' ^PoaCp.  ' 

L 3a)Po,3c,a''  J 

1 ' 2po,Cp3 ) 

rCwC,(Ya-1)] 

/ (DK, 

\ 

1 topo^ab^  J 

V 2Po,Cp, 

1 

+ I 


^C',K,(Ya-1) 

6po«Po,^c,Cp,a2b2 


+ 1 


r 11^^c,(Y.-1)  -1  / a)K, 

L 3(jopo^Po,2c,a3b2  J V 2po„c„,  / \ b3  / 


2POaCpa 


r 2^C^C,(Ya-1)1  / («JK, 

I w^PowPota^b^  J V 2po,Cp,  ^ V b3 


J £li25C«K,(Ya-1)l  r,  , 
I 2Po,3C,3Cppa3  J I 


f a)^2c„(Ya-1)l 

r+i 

r ^C«(Ya-1)l 

/ WKa  /'ll  , 

1 po,3c,3a3  J 

UpOaCp,  ) 

L 6Po,3C,b3  J 

V2po,Cp3  j 1 ^ a3  j 

(B-58) 

The  terms  in  equations  B-57  and  B-58  have  been  arranged  so  that  the  real  and  imaginary  inertial,  viscous, 
thermal,  and  viscothermal  terms  are  grouped.  These  terms  will  be  designed  as  the  m-,  q-.  ^q-,  im-,  i^-, 

iq-,  and  i^q-terms,  respectively.  In  equation  B-57,  the  first  five  terms  are  m-terms.  These  are  followed  by  two 
^-terms,  five  q-terms,  six  ^q-terms,  four  i^-terms,  five  iq-terms,  and  six  i^q-terms.  in  equation  B-58.  the  first 
four  terms  are  ^-terms.  These  are  followed  by  five  q-terms,  six  ^q-terms,  five  im-terms,  two  i^-terms,  five  iq- 
terms,  and  six  i^q-terms. 

The  next  step  is  to  simplify  each  group  of  terms  for  Y and  (Y-Z).  However,  examinations  of  equations  B-57 
and  B-58  indicates  that  the  powers  of  the  parameters  in  the  Y and  Z terms  are  never  identical.  Thus,  the  terms 
in  Z can  be  examined  independently  of  Y,  since  the  subtraction  of  Z from  Y can  not  result  in  the  reduction  in  the 
order  of  magnitude  of  any  terms.  After  each  group  of  terms  is  compared  and  simplified,  the  expressions  for  Y 
and  Z will  be  further  simplified  by  comparing  all  remaining  terms.  Then  the  subtraction  (Y  - Z)  will  be  per- 
formed and  that  expression  simplified. 

Comparison  of  the  terms  within  the  groups  in  equation  B-57  yields: 


/ (DC, a \ 

/ (D3^3b  \ r llPo,c,3  1 

^ - 1 1 

V 3pp,b3  ) 

V3Po,3C,3a3/  1 9Po,,C,3  ' 

r b3 ’ J 
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' » 


r 2^c».c,(Ya  - 1) 

1 ( r 

f 1l0)^  \ 

1 o)2po^Po,a3b2 

J \ 2Po3Cpp  / 

6Po,q2  ) 

3PoaCaL\  /'  1 ^ 


3p0aC,"a 


)-,] 


(- 


& (-S^) 


SpowPc^Qa^b^ 

_ : r 2^C^Cl(Ya-1)  1 / COKa  \ '"  / . _ HtD^  \ 

L c»"Po»Po,a^b"  J \ 2poaCp3  / ( 6po,c,2j  (B-61) 

Then  a comparison  of  all  terms  in  equation  B-61  yields: 

/ PcPoaCa"  \j  ^ \ r C^Po,(Ya-1)-[/  ^^a  ^ _£f_\ 

\ C(  / V 3a)Po^a^b2  / I copo^ab^  jV2po.Ca.  / V’  Po,(i)2a3/ 


3a)Po.a2b2 

3p0aCa= 


'(-3?^)  [ (■po,w"a2 

. , (■  

'-9Po.Po,c,2a2b2 


)(- 


copo^abJ  JV2poaCp3 

3PoaCa"a  / J 


)-[ 


. f ^wPo((Ya ^ ) 


topo^ab^ 


i(^)  ( 


- 2s 
1 + 


Po,w2a3 


■)  • (B-62) 


Examination  of  equations  B-60  and  B-62  shows  that  all  the  viscothermal  terms  are  negligible. 
Subtraction  of  equation  B-62  from  B-60  and  a subsequent  comparison  of  terms  yields: 


PcHPoa^-a^ 


3Po,Cw(Ya-1)-|/  /,  , 3s  \ 

JUpOaC^/  V Po,W"a3/ 


0)2  a2 


[•  -^PoiC 

L Pow 

* '(^)  l(^) ( ’ -ro)  - - (ra-)] 

Equation  III-32  indicates  that  I B»  is  required  to  determine  o.  This  implies,  frorr  equation  B-48,  that  I Y I2 
is  required.  If 

Y = Y,  + iY,  , 


then 


|Y|2  = Y,2  + Y2  , 


(B-64) 

(B-65) 


where  Y,  and  Y,  are  the  real  and  imaginary  parts  of  Y,  respectively.  Thus,  in  light  of  equation  B-65,  it  is  possible 
to  simplify  equation  B-60  by  comparing  magnitudes  of  real  and  imaginary  terms.  This  comparison  yields: 


(B-66) 
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the  swimbladders  of  small  fish.  A new  model  of  a swimbladder-bearing  fish  has 
been  developed  in  order  to  provide  improved  predictions  of  the  resonant 
frequency  and  acoustic  cross  section  of  such  a fish.  Development  of  a new 
model  was  undertaken  because  comparisons  of  the  predictions  of  previous  models 
with  experimental  data  show  these  models  to  be  inadequate.  Primarily,  the 
experimental  data  indicate  that  deutiping  in  fish  tissue  is  consistently  greater 
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than  that  predicted  by  previous  models.  In  addition,  in  several  instances, 
resonant  frequencies  have  been  measured  which  are  significantly  higher  than 
can  be  accounted  for  by  these  models. 

The  new  model  consists  of  a small,  spherical  shell,  enclosing  an  air  cavity, 
water.  The  shell  is  a viscous,  heat-conducting  Newtonian  fluid,  with  the 
physical  properties  of  fish  flesh.  The  interface  between  the  shell  and  the 
cavity  supports  a surface  tension.  The  shell  is  insonified  by  a harmonic 
plane  compressional  wave,  whose  wave-length  is  large  compared  to  the  shell 
diameter. 

The  solution  to  the  problem  consists  of  obtaining  the  amplitude  of  the 
scattered  compressional  wave.  This  was  accomplished  by  first  determining 
the  appropriate  wave  equations  from  the  linearized  equations  of  motion. 
Eigenfuction  solutions  to  these  equations  were  then  obtained.  However, 
theoretical  evidence  indicates  that  only  the  fundamental  mode  is  an  important 
contributor  to  volume  reverberation,  so  that  higher  modes  were  neglected. 
Application  of  the  boundary  conditions  then  led  to  the  required  solution  for 
the  amplitude  of  the  scattered  compressional  wave.  Due  to  the  complexity  of 
this  solution,  order  of  magnitude  ainalyses  of  the  various  terms  were  con- 
ducted in  order  to  obtain  a simplified  expression.  Equations  for  both  the 
resonant  frequency  and  acoustic  cross  section  were  then  determined  from  this 
simplified  expression. 

A comparison  of  the  results  of  the  new  model  with  experimental  data  indicates 
that  the  new  model  constitutes  a definite  improvement  over  previous  models. 
The  new  model  can  predict  the  high  values  of  damping  and  elevated  resonamt 
frequencies  that  previous  models  could  not.  The  comparison  indicates  that 
the  model  is  most  accurate  for  fish  in  which  tension  in  the  swimbladder  wall 
has  a minor  effect  on  resonant  scattering.  This  includes  the  fish  which  are 
of  interest  in  studies  of  volume  reverberation  and  therefore,  the  new  model 
will  be  of  considerable  value  in  such  studies. 


